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Summary. Pulsars are the quintessential dogs that don't bark in the night - their 
observed loss of rotational energy mostly disappears into the surrounding world 
while leaving few traces of that energy loss in observable photon emission. They 
are the prime example of compact objects which clearly lose their energy through a 
large scale Poynting flux. 

I survey recent successes in the application of relativistic MHD and force-free 
electrodynamics to the modeling of the pulsars' rotational energy loss mechanism 
as well as to the structure and emission characteristics of Pulsar Wind Nebulae. 
I suggest that unsteady reconnection in the current sheet separating the closed 
from the open zones of the magnetosphere is responsible for the torque fluctuations 
observed in some pulsars, as well as for departures of the braking index from the 
canonical value of 3. I also discuss the theory of high energy pulsed emission from 
these neutron stars, emphasizing the signiflcance of the boundary layer between the 
closed and open zones as the active site in the outer magnetopshere. I elaborate 
on the conflict between the models currently in use to interpret the gamma ray 
and X-ray pulses from these systems with the electric current flows found in the 
spin down models. Because the polar cap "gap" is the essential component in the 
supply of plasma to pulsar magnetospheres and to pulsar wind nebulae, I emphasize 
the importance of high sensitivity gamma ray observations of pulsars with core 
components of radio emission and high magnetospheric voltage, since observations 
of these stars will look directly into the polar plasma production region and will 
probe the basic plasma parameters of these systems. I also discuss the current state 
of understanding and problems in the shock conversion of flow energy into the spectra 
of the synchrotron emitting particles in the Nebulae, emphasizing the possible role 
of heavy ions in these processes. I comment on the prospects for future developments 
and improvements in all these areas. 



1 Introduction 

This paper was originally entitled "Pulsar Emission: Where to Go" in the 
theory of pulsar radio and high energy emission. "Where to Go" on these 
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topics depends greatly on where one thinks the energy and the particles are. 
Therefore, I mostly discuss questions of dynamics, and return to the emission 
physics through the lens of the successes and problems of dynamical models 
of rotation powered pulsars' magnetospheres. Also, I emphasize results and 
problems common to all pulsars. I will give an impressionistic rather than 
a comprehensive review, more in the spirit of setting goals as I see them 
rather than providing a scholarly survejQ. I have focused on issues that can 
be addressed by timing and by high energy photon observations. For a survey 
of the current status of models of pulsar radio emission, see the article by 
Usov in these proceedings. 

2 Pulsar Electrodynamics: Follow the Energy 

Astrophysical understanding comes from using observation and theory to find 
and follow the flow of energy, mass and momentum in the macroscopic sys- 
tems of interest. For pulsars, this has been a challenge. They are quintessential 
"dogs that don't bark in the night" , with only a small fraction of the energy 
they broadcast into the Universe appearing in directly observable forms. The 
interpretation of the regular pulse periods from sources distant enough to 
require stellar and sub-stellar luminosity in the radio (and in the infrared, 
optical. X-ray and gamma ray) discovered by the radio astronomers led im- 
mediately to the understanding that the observed periods are the consequence 
of rotation of massive stellar flywheels (neutron stars) . The steady lengthen- 
ing of the pulse period, shown in Figure (TJ led immediately to a permanently 
successful model of that spindown, the electromagnetic torques exerted if the 
stars are sufficiently well magnetized. 

One can readily estimate the magnitude of such torques from the observa- 
tion that rotation of a stellar magnetic field B induces a poloidal electric field 
of magnitude E ~ {ilr/c)Bp, with /? = 2tt/P and Bp the poloidal magnetic 
field - from the point of view of the torque, that field is well approximated by a 
dipole with dipole moment ^ = R^Bp. The winding up of the magnetic field as 
the conducting star rotates requires the existence of a toroidal magnetic field 
of magnitude B^f, ^ {Qr/c)Bp . This E field corresponds to energy loss in a 
Poynting flux cE x B/Att. If the electromagnetic energy density exceeds all the 
material energy densities, one obtains the total energy loss Er and therefore 
the torque J — En/fi by summing the Poynting flux over a sphere of ra- 
dius expected to be comparable to the light cylinder distance Rl — c/H, 
where the electromagnetic inertia B^/Attc^ causes the poloidal field to depart 

^ Much of what I have to say derives from collaborations, most recently with Elena 
Amato, Phil Chang, Niccolo Bucciantini, Eliot Quataert, Todd Thompson and 
especially Anatoly Spitkovsky; in earlier years, with Ted Scharleman, Bill Fawley, 
Colin Norman, David Alsop, Don Backer, Brian Gaensler, Yves Gallant, Vicky 
Kaspi, Bruce Langdon, Claire Max and Marco Tavani. However, I am solely re- 
sponsible for the views expressed in the subsequent pages. 
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from the imposed stellar (dipole) field by an amount on the order of Bp itself. 
Then i?^ is the smallest radius where i?^ becomes comparable to Bp, and 
Bp w IJ-/R\- Therefore 



47r \ J \Ra 

For radio emitting Rotation Powered Pulsars (RPPs), stars are known with 
Eu — IQfl from as small as 10^° ergs/s to as large as 10'^^ erg/s - / is the 
stellar moment of inertia, / w 10^^ cgs for currently acceptable equations of 
state for neutron stars. 
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Fig. 1. Observed RPP periods and period derivatives, from 
iKaspi. Roberts fc Harding (|2006l ) . "X" marks a pulsar with P, P measured 
from X-rays as well as radio observatons. The line <P = 10 Volts = 1 TV is the 
locus in the PP diagram where the rotation induced voltage drops to 10^^ V, 
clearly marking a boundary beyond which pulsar emission is unlikley 



Modeling of RPPs has one great advantage over modeling of other compact 
objects - observations of P, P determine the energy supply, to within the 



4 Jonathan Arons 



uncertainties in the moment of inertia. In contrast, modeling of accreting 
black holes always suffers from major uncertainty as to whether the systems 
are, or are not, accreting at a well determined rate, e.g., the Eddington limit. 
This fact runs through much of what I discuss below - with the total energy 
budget known, the effort turns to aspects of the machine's physics at a level 
of sophistication not sustainable in many other aspects of compact object 
physics. 



2.1 Force-Free Model: Heuristics 

Expression ^ makes no reference at all to charged particles, and indeed the 
first theories of RPP spindown (some invented before RPPs were discovered) 
invoked the electrodynamics of a vacuum rotator as an explan ation of the 
observed P ( Deutsch 1955t Pacini 1967 : Ostriker &: Gunn 1969f ). Except for 



geometric factors, vacuum theories yield expression ([T]), but with the special 
addition that as the angle i = Z(/x, f2) becomes small, so does the torque, 
in proportion to sin^ i. Application of this model to the spindown data for 
normal RPPs yields dipole moments on the order of lO'^" cgs, corresponding 
to surface dipole fields = fiR^ ^ 10^^ Gauss for "normal" neutron star^j 
Vacuum models have large electric fields parallel to the magnetic field at 
the stellar surface, a fact which led Deutsch (in the context of the oblique 
rotator z 0) to suggest that a vacuum rotator has to form a charged mag- 
netosphere, as charged particles move from the surface to short out E ■ B. 
Simultane ously with the appe arance of the vacuum torque models after RPP 
discovery, Goldreich fc Juliai] f|1969.,) independently made the same observa- 



tion in the context of the aligned rotator. They went further to suggest that 
a charge separated outflow forms, creating a conduction current J = rjnv, 
where ry^j is the charge density required to force = E ■ B / B to zero, 
-qji — ~fl- B/2tic + relativistic corrections. They also introduced the idea 
that the magnetosphere is "force-free", that is, the electromagnetic energy 
density is so large that all inertial, pressure and dissipative forces can be ne- 
glected, a concept consistent with the fact that RPPs are non-barking dogs - 
the large energy loss manifested in spindown does not appear in any radiative 
emission associated with the magnetosphere (here defined as the region inte- 
rior to Ra, probably k, Rl-) Conceived of as a system which is strictly steady 
in the corotating frame - after all, pulsars form superb clocks, therefore the 
rotating lighthouse picture should apply, which it does, at least to averages of 
many pulses - the flow of the charges decompose into any velocity parallel to B 



"Normal" means neutron stars discovered via their "normal" radio emission, an 
obvious selection effect. More recent discoveries, of millisecond radio pulsars and 
of X-ray selected objects, have revealed neutron stars with magnetic moments 
from ~ 10"^'^ CR S down to "zero" in the X-ray burs t sources, which effectively means 



10^® cgs (|Kaspi. Roberts fc Hardi"ng |2006| ). In particular. X-ray cyclotron 



^l < 

lines confirm the existence of 10^^ Gauss surface fields. 
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plus rotation f2 x r. The same charge density and velocity dec omposit i on ap - 
ply to the magnetohydrodynamic (MHD) model introduced bv lMichell ( 19691 ) 
in the same year, with the difference that the MHD model assumes a density 
large compared to rjR/q. Charge separated/MHD outflow and magnetospheric 
dynamics occurs for densities equal to or less than/greater than 
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if the particles have Larmor radii small compared to r. Here Pioo = P/100 
msec and /iao = /i/10'^'^ cgs. Interpreted as a particle outflow, this density 
corresponds to a particle loss rate 



Nr. 



. n2 \Vr{Rl 



— = — = 2.7 X 10^2^ s" 



(3) 



The poloidal electric current In is the current expected such that the induced 
magnetic field becomes comparable to the dipole field at the light cylinder. 
In Goldreich and Julian's charge separated picture of the aligned rotator, 
the charges in In, composed of the charges in the fully charge separated 
plasma flowing parallel to the poloidal magnetic field, provides the support 

for i?0 and the Poynting fiux. <P = \J Er/c = 1.3 x 10^'''/X3o/Pioo Volts is the 
magnetospheric electric potential. 

Observati ons of the synchrotron emission from young pulsar wind nebulae 
CPW Ne'l ffi.o.lRees fc Guiinl ll974l: iKennel fc Coronitilll984l : ICxallant fc Aro^ 
1994 iGaensler et aLll2002[ ) reveal particle injection rates N (in the form of 
electrons and - probably - positrons) corresponding to densities in a wind 
outflow n = N /ATir^c (at distances much larger than Rl) a factor of lO'^ 
and more larger than the density of the charge separated flow predicted by 
the charge separated wind idea. Thus the nebular observations suggest MHD 
models with a quasi-neutral plasma (which can only be electron-positron pairs, 
see ^J3]), appear to be a good starting place for understanding these systems. 



2.2 Force-Free Model: Results 

Thus the simplest idea is that a dense plasma exists everywhere in the mag- 
netosphere and beyond, with the plasma energy density much lower than 
B'^/Stt - for the young, high voltage pulsars, plasmas with energy density 
remotely comparable to that of the EM fields and still under the rotational 
control of the stars would lead to pulsed photon emission orders of magnitude 
greater than what is observed. The force free idea was elegantl y formulated in 
the early 70s in the "pulsar equation" for the aligned rotator (|Miche]|ll973al : 
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Scharlemann fc Wagoneilll973f ). a variation of the Grad-Shafrono v equation 



familiar from the theory of magnetic confinement (|Batemanlll978l ) 



d^^ a^VA / 1 , .dl 



Here ip is the magnetic flux, with the poloidal magnetic field related to ip by 
Bp = —w~^4' X Vi/i, vj is the cylindrical distance from the rotation (z) axis, 
while the toroidal magnetic field is — I{ijj)/'cu, with / the current enclosed 
within a circle around the z axis of radius zu. These fields are supported 
by charge and current densities all derivable from once a solution of Q 
is determined. The electromagnetic structure of this non-pulsing model was 
(and is) thought to capture the essence of what is needed to make a full, 
oblique rotator model. The model says nothing about particle energetics - 
thus it, and its oblique rotator descendants, provides a geometric platform 
and an accounting of the dominant electromagnetic energy flow tapped by 
the subdominant dynamical processes that lead to the observable emissions. 
In particular, it does provide a basic model for the invisible processes that 
lead to pulsar spindown. 

Solution of (HI) in the simplest relevant case (a star centered dipole with 
rotation axis parallel to the magnetic moment) has taken a remarkably long 
time. Solutions appeared immediately for a) a strictly co-rotating magneto- 
sphere I{tp) = 0, which is not relevant since it does not spin down (this dog 
really doesn't bark!) and implies particle motions faster than the speed of 
light at w > Rl a n d b) a star centered monopole, with an elegant result 
obtained by iMichell (|l973d) whose most important element is the poloidal 



current function 

^<«^-|(^)' 

where ipo « '!rw'^^p{2^/ Rl) « Rl){Rl/Ry) is the open magnetic flux in 
one hemisphere of the monopole - Ry is the equatorial radius of the Y point 
in the magnetic field which marks the largest extent of the closed magnetic 
field hues in the rotational equator, Wcap ~ {R^,/ RyY^'^ is the cylindrical 
radius of the magnetic polar cap, and Vo = Ry'^^ — y^jRy- Finding these 
solutions required inspired guessing of /('(/')■ In the years between 1973 and 
1999, many attempts were made to solve ([4]) by guessing various forms for 
/(■0) and applying ever more clever analytic techniques t o this fundame ntally 
non- linear model. None yielded anything credible - see Mestel ( 19991 ) for a 
summary of much of this work. 



The situation changed when IContopoulos et al. ( 19991 ) took seriously the 



nonlinear eigenvalue and eigenfunction character of (j4]) and its associated 
boundary and regularity conditions and successfully applied an iterative nu- 
merical technique to find B and to produce a result with E ■ B = 
and — < everywhere - the latter condition is required if the model 
is to be taken seriously as a representation of a physical magnetosphere. 
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since the E x B velocity must be subluminal, for a physical modeH. This 
solution, in which the last closed field line is assumed to have equatorial 
radius equal to Rj,, ha s been repr o duced with incr easing numerica l accu- 



racy b v lGruzinovl (|2005f) . lMcKinnevl (|2006l ) . iTimokhin (2006) and Spi tkovskv 
( 2006f ). It exhibits a number of long expected features (jMichellll975( ). In par- 
ticular, the last closed field line has a Y-type neutral point on the equator, 
with return current flowing (mostly) in an unresolved current sheet along 
the boundary of the closed zone, then extend ing as an equato rial current 
sheet to radii w > Ry- Figure [21 taken from Timokhin ( 20061 ). shows the 
poloidal magnetic geometry of the ahgned rotator. The solution illustrated 
has '0n = 1.23u,/ Rr,, in th e case Ry/Rt, = 0.992, in excellent agreement with 



r-^rA'- I 

Gruzinovl (j2005l ) and with IContopoulos et al\ (|1999l ). who assumed Ry — Rl 



exactly. All authors agree on the spindown energy losses of the aligned rotator, 

Er = kQ^n'^/c^, k = l±p -l. 

As predicted bv lMichel ( 1974 ). the asymptotic structure (r — \Jv3^ + ^ 
Rl) approaches that of the (split) monopole, as appears most clearly in 
McKinney's and Timokhin's results. Thus the poloidal current flow is almost 
that of the monopole , a point discussed furth er below. Also, a s shown by 
Goodwin et^ (|2004[) . IContopoulosI (|2005l ) and iTimokhh] (|2006l ). the steady 
state force free magnetopshere has a whole range of possible solutions, param- 
eterized by Ry /Rl < 1- 




Fig. 2. Left: Cartoon of the aligned rotator's magnetosphere, showing the primary 
polar current and the return current flowing along the separatrix a the equatorial 
current sheet. The anti-aligned case, with /i anti-parallel to fl is shown. The aligned 
case has the same topology, with the sign of the current flows reversed. Right: 
Field lines (magnetic flux surfaces) of the aligned rotator solution, for the case 
Ry = 0.992Rl 



^ — > is possible, in principle. However, in the absence of losses particles 
then accelerate to energy ~ in distances not greater than Rl In the younger 
pulsars, the acceleration becomes radiation reaction limited, implying radiation 
emission from RPPs far in excess of what is observed. 
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By solving the time dependent force free equations, ISpitkovskv (I2OO6I) 



showed that the force-free magnetosphere evolves with Ry — ^ Rl, starting 
from a static vacuum magnetic dipole on a star instantaneously set into rota- 
tion with angular velocity f2; at t = 0, the electric field on the stellar surface 
was set equal to —{f2 x r) x B. The rate of approach of Ry to R^ depends 
on the artifi c ial res istivi ty used to control s ingula r behavior at current sheets. 
Komissarov ([2006!) and Bucciantini et al. ( 2006^ found similar results using 



a relativistic MHD model {i.e. inertial forces included). Spitkovsky's method 
allowed the current sheet to have an arbitrary shape. Thus, he also succeeded 
in finding the force free solution for arbitary i; the resulting 3D model of the 
magnetic field appears in Figure[3l Within the assumption of a magnetosphere 
everywhere filled with plasma of density sufficient to short out parallel elec- 
tric fields (and no physics that might support such electric fields in a plasma 
of density greater than \r]n\/e), a full solution for the electromagnetic struc- 
ture of the RPP's magnetosphere (both aligned and oblique rotators) is now 
available, after 38 years of discussion. Spitkovsky's result. 



Er = kf-^(l + sin^ i), k=l±0.1, (6) 

contains all the aligned rotator studies as a special case. As is clear from Figure 
[31 the magnetic topology of the oblique rotator (closed field lines terminating 
at a Y line, current sheet extending from the Y line separating regions of 
oppositely directed field in the wind) is a rotationally distorted version of the 
simpler aligned rotator geometry. 




Fig. 3. Left: A snapshot of a force free simulation of a RPPs magnetosphere, for 
r < 2Rl (from Spitkovsky 2006). Right: Total current (c/47r)V x B, the sum of 
conduction and displacement currents. 
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2.3 Beyond the Force-Free Model: Plasma Sighs and Whispers 



Until recently, the magnetosphere was assumed to have Ry — Rl, an assump- 
tion consistent with the observation that radio emission from low altitude 
appears to occupy a polar flux tube which, if modeled as being in a static vac- 
uum dipole's geomet ry, is bounded by a closed zone w hich appears to extend 
to the light cylinder ( Rankinlll990t Kramer et al. 1998 ). i.e., have a polar cap 
opening angle 6^ — {R^/RyY^'^ with Ry = Rl', of course, the observations 
and the simplified model do not come close to proving that Ry = R^ , or even 
that static dipole model for the B field is accurate all the way to r = Ry. But 
if Ry / Rl < 1 is some constant, then since Rl increases as the pulsar spins 
down, there must be net conversion of open field lines to closed field lines 
(transfer of open magnetic flux to closed flux) on the spindown time scale. 
This topological change requires rcconncction and a violation of ideal MHD, 
at least in local regions - the likely culprit is the Y -line and the current sheet, 
as has been observed bv IContopoulod ( 20051) and Contopoulos fc Spitkovskv 
( 2006|) . Having made this observation. I Contopoulod ( 2005f ) suggests that for 
pulsars near the death line in Figure [H reconnection proceeds sporadically - 
the magnetosphere "coughs" - because of failure in the supply of plasma from 
the polar cap accelerator and pair creation region - see ^ He applies this idea 
to a scenario for the m ajor outburst observed in the magnet ar SGR1806-20 
on December 27, 2004 (|Hurlev et aLlbOOsHPahner et aZ.ll2005[) . 

In fact reconnection is likely to have an unsteady, "bursty" character for 
all pulsar^ - the magnetosphere should be noisy at some level all the time. 
Figure [4] shows a snapshot of the m agnetosphere of a rotati ng neutron star 
with large plasma supply, taken from Bucci antini et al. ( 20061 ). This relativis- 
tic MHD (not force-free) model was designed to represent the wind from a 
young, rapidly rotating magnetar, with the wind driven by the enormous 
thermal pressure at the neutron star's surface. The wind, formed by plasma 
flowing out on open field lines, converges on the equatorial current sheet. 
That convergence causes driven reconnection, operating in a bursting mode - 
the sheet forms "plasmoids" islands of reconnected poloidal field with closed, 
0-point magnetic topologjO which flow away at the local Alfven speed, ~ c. 
Numerical resistivity in the code provided the dissipation required to allow the 
transformatio ns of field topo logy shown - the GEM study of non-relativistic 
reconnection (jBirn al\ 120011 ) demonstrated that any non-ideal effect allows 
rapid driven reconnection with inflow velocity into the separatrix (the current 
sheet separating the closed field from the open field regions in Figures [2] and 
[3|) being (0.1 — 0.2)va, outside the restrictive bounds of incompressible MHD 
with uniform resistivity. Note that the reconnection sporadically transforms 
the Y-line into a X-line, with the current she et then cont aining a series of 
dynamical X-lines, all leaving the star - ,Lvubarskv 3 (|l990h objection to the 



* Reconnection measured in the laboratory and in space plasmas, and observed in 

solar plasmas, does occur with bursty, often explosive, behavior. 
^ In these axisymmetric models, the islands are magnetic torii. 
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formation of a stationary X-line (the field lines inside the separatrix on the 
open side of the X-line are not anchored to the magnetosphere) is answered 
simply by the fact that the plasmoids indeed are not anchored to the magne- 
tosphere and fly away, but constantly reform. The Poynting flux was found to 
be time dependent, fluctuating around the mean by ~ 30%. 

The Bucciantini et at model was not designed specifically for classical 
RPPs, or for magnetars in their currently observed state - future work on 
reconnection in configurations with electron-positron plasma, where the Hall 
effect, so important in the nonrelativistic studies, is absent (and is replaced ei- 
ther by pressure anisotropy as th e facilitator of rapid reconnection, as pointed 
out bv lBessho &: Bhattachariegi2005l . or by particle inertia), in a state sug- 
gested by the models of plasma supply in RPPs' magnetospheres, are required 
to quantify this "noisy magnetosphere" picture. In particular, extending such 
modeling to the full 3D rotator in strongly magnetized MHD has not yet been 
done. Also, it remains to be demonstrated that noise in the current sheet at 




Dip: magnetic field Dip: magnetic field + J 




Fig. 4. Left: Magnetic structure of a relativistic magnetosphere with mass outflow, 
in a case when the magnetic pressure at the light cylinder exceeds the relativistic 
plasma inertia by a factor approaching 20, a record high for MHD calculations. 
The contours represent poloidal magnetic field lines, while colors represent the ratio 
B^/Br- Reconnection have been artificially suppressed. Middle: Snapshot of the 
magnetic structure when the reconnection components of the electric field are not 
suppressed. Reconnection occurs because of numerical resistivity introduced by the 
finite difference scheme. Plasmoids emerge along the equatorial current sheet, grow 
and flow out at the local Alfven speed (~ c). Right: Blow up of the plasmoid 
structure. 
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and beyond Ry communicates back to the inner magnetopshere and the star, 
through (kinetic) Alfven waves traveUng back along field lines at and near the 
separatrix (see 

But if this picture does apply to RPPs, it has a number of consequences 
for observables and outstanding questions, some of which I touch on further 
in J|3l From the point of view of the basic energetics embodied in spindown, 
the fluctuating Poynting flux may imply a fluctuating torque. Noise in pulsar 
spindown has been known since the early days - it limits the ability to time 
pulsars coherently. If the magnetic field interior to but near Ry fluctuates 
by tens of per cent on time scales comparable to the rotation period, and 
these fluctuations represent variations in the poloidal current that communi- 
cates stress to the star, then the torque is noisy with magnitude the same as 
is inferr ed from representing the observed random walks in the rot ation fre- 
quency (jHelfand et al\ Il980l: ICorded Il980t ICordes fc Helfandl Il980l) as bein g 
the consequence of white noise in the electromagnetic torque ( Aronsl[l981b - 
reconnection may provide the mechanism for magnetopsheric instability and 
torque fluctuations that was not specified in this early attempt at scenario 
building) . 

Magnetos pheric noise opens the p ossibility that Ry/Ri, evolves. One can 
readily show (jBucciantini et aZ.ll20Q6l ) that the braking index, defined as n = 
flfljfP' (with averages indicating the usual average over subpulses taken in 
measuring pulsars' periods, here taken to be the same as an average over 
plasmoid emission and torque fiuctuations), in a magnetosphere with evolving 
Ry / Rl but fixed /i and i, is 



nn 



= 3-^2 



a In [\ 



Ry 
Rl 



51n/2 



(7) 



Braking indices less than 3 (jLivingstone et a/.ll2005l and references therein) 
thus may indicate a progressive lag of the closed zone's expansion, measured 
by the radius of the Y line, behind the expansion of the light cylinder as a 
pulsar spins down. 

This is hardly the only thinkable explanation of n < 3. Magnetic moment 
evolution has long been advocated as the o rigin of small braking indices, 
going back to the crustal field growth model of B landford et al. (1983) - which 
doesn't actually work in those authors' formulation, the threshold for growth 
set by crustal resistivity is too high - to Ruderman's model for growth of 
M± = M siu i oc ^/P due to interaction of interior ma gnetic flux tubes with the 
quantized vorticity of the superfluid interior (e.g., Rudermanl 20061 ). When 
combined with expression ([S]), this model yields 



3 4^ tan^ io 



1 



2^tan2 io' 



(8) 



where (Iq, io are the angular velocity and obliquity at the time when the 
arrays of quantized magnetic flux tubes and vortex tubes have both formed. 
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thought to be perhaps lO'^ years after the neutron star's bu'th. This model 
can produce any braking index between 2 and 3. 



Pulse energy Single— pulse seguence 



5000 
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Fig. 5. A series of individual pulses from PSR B0943+10 (center). The average pulse 
is shown in the bottom panel, as a function of pulse number and longitude, with 
360° or longitude corresponsing to one rotation of the s tar. This star shows or g anize d 
drifting of the pulses through the pulse window. From lOeshpande fc RankinI (|l999l ). 



Fluctuations and oscillations in the corotating frame of the currents at 
the light cylinder also offer a possible explanation for the long known fact 
that pulsars are flickering lighthouses. The well-known pulse stability that 
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allows exquisitely precise timing applies to average pulses, formed by sum- 
ming hundreds to thousands of individual pulses. However, individual pulses 
arrive with at varying times within the pulse window, usually at random in 
those pulsars with "core" emission characteristics, and either at random or 
with an organized drift of the arrival times t hrough the p ulse window, in 
stars with "conal" emission characteristics (see Rankin 19831 for pulsar beam 



classification) . Figure [5] shows an example of pulse to pulse variability. Typ- 
ically one or at most two pulse components are within the pulse window at 
any one time, suggesting the individual pulse variability time is on the order 
of the rotation period. That time scale is consonant with the Alfven wave 
transit time from the low altitude emission region to the radius of the Y-line, 
where current fluctuations are formed, a coincidence suggesting the subpulse 
variability (both random and drifting) is a con sequence of c urrent variations 
created by activity in the outer magnetosphere ( Arond 1981b ). In this picture, 



drifting and chaotic subpulses are both the consequence of the same dissipa- 
tive dynamics of the currents coupling the magnetosphere to the wind, with 
drifting subpulses reflecting limit cycle behavior of magnetospheric reconnec- 
tion while chaotic subpulses represent a more random, bursty behavior of the 
field lines topological changes. Objects such as PSR B0943-I-10, in which tran- 
sitions fro m organized drifting to chaot ic single pulse behavior and back are 
observed ( Rankin fc Sulevmanoval [20061 ). are particularly teUing laboratories 



for investigation of the connection between current flow and pulsar emission, 
and thus offer insight into magnetospheric dynamics. 

2.4 Electrospheres? 

There is, however, a ghost hiding inside the force-free/MHD magnetospheric 
machine. These theories assume that a plasma dense enough to enforce £'|| = 
is present everywhere in the magnetosphere, an assumption which relies upon 
the pair creation physics summarized in Jj3l Pair creation assumes relativistic 
beams contribute a substantial fraction - possibly all - of the electric cur- 
rents embodied in the force-free and MHD models, since only these beams 
can (plausibly) emit the gamma rays that convert to e^. One can think of 
pair creation as an instability of the current flow originaUy hy pothesized in 
the fully charge separated scenario of lGoldreich fc .Tulianl (|l969h . But, as was 



recognized not long after the charge separated outflow scenario was suggested, 
the charge separated wind model must fail, so long as charged particle flow 
across field lines is forbidden - many field lines of a dipole (not a monopole) 
must pass through a surface where f2 ■ B = 0. The charge (and plasma) den- 
sity of the charge separated medium on the exterior (larger) radius side of this 
"null surface" has sign opposite to that of the plasma that can be supplied 
from the stellar surface by particle motion parallel to B. The plasma in this 
exterior region has no source, if the only allowable charged particl e motions 



are sliding along the magnetic field plus bulk flow E x B drift ( HoUowavl 
197.1 11975^ . Thus, one expects such a magnetosphere to open large gaps, and 
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probably have no charged particle wind - certainly no wind with particle flux 

greatly in excess of TV^ = c<^/e. 

Such "electrospheres" (IKrause-PolstorfF fc Michellll985r ) do not appear to 



be unstable to pair creation (jPetri et aLll2002a[ ). thus do not collapse to the 
plasma filled state hypothesized in the force free models, by pair creation 
alone. However, they are unstable to cross field transport even without pair 
creation. A large gap separates the equatorial regions of the electrosphere from 
the stellar surface, leading to differential rotation of the equatorial plasma. 
This differential rotation is subject to the di ocotron instability - a vari ation o f 
a Kelvin- Hclmholtz instability (Spitkov skv fc Arondl20 02: Pctri__ej_Qi. 2002bl ). 



Simulations (jSpitkovskv fc Aronsll2002l : IPetri g/Hiooal : |Petrill2007h suggest 
the resulting time dependent (in the corotating frame) E x B drifts create 
cross-field "diffusion" which may relax the charge separated magnetosphere 
to something approximating the state envisaged by Goldreich and Julian. 
Expansion of the equatorial plasma is illustrated in Figure [5] 

Bottoms-up models based on these results have not been investigated. 
The fact that young pulsars supply their nebulae with particle fluxes greatly 
in excess of Nr lends support to the perhaps more practical view that the 
filled magnetosphere model has consequences in reasonable accord with ob- 
servations of high energy pulsed emission and of pulsar wind nebulae, thus 
deserves the main focus of scientific attention. Such models may have charge 
separated current flows possibly unstable to pair creation. 

The force free solutions and their possible extensions have a number of im- 
plications for emission models, and for the pair creation models that underlie 
the emission physics. 



2.5 Magnetic Geometry of Radiating Layers 

The oblique rotator solution determines a polar cap/polar flux tube size and 
shape. These are noncircular, and have centroid displaced from the magnetic 
axis (Spitkovsky, personal communication). Such changes in magnetic geom- 
etry from the conventional assumptions need to be folded into radio beaming 
and polarization, and polar cap X-ray emission models, which often invoke ad 
hoc changes in polar cap size and s hape, perhaps created by surface anomalies 
in the surface magnetic field {e.g., Melikidze et al. 2006I ). in order to explain 



departures from the simplest, static dipole geometry. The need for such extra 
parameters in the models needs to be evaluated in the context of realistic 
magnetic geometry of the rotating dipole. 

The surface of last closed field lines (the separatrix) and of the return 
current flowing along that separatrix has been determined within the force- 
free approximation. Particle acceleration in gaps (regions of low density where 
a parallel electric field = E ■ B/B forms because of charge starvation 
below the Goldreich- Julian density) on the open field lines close to this surface 
has been advanced as the o rigin of the pulsed gamma rays observed from 
a small number of pulsars (jMuslimov fc Hardin g| 120041 : iHirotanil lioiil and 
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references therein) by the EGRET experiment ( Thompson 20001 ) and by other 
high energy detectors, with a substantial increase in t he population observed 
expec ted with the launch of the GLAST telescope ( McLaughlin fc Corded 
20001) . 
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Fig. 6. Spreading of the equatorial charged disk in the electrosphere of an aligned 
rotator un der the influence of the dioc otron instability, from a PIC simulation of the 
flow, from lSpitkovskv fc Aror3 (|2002l ). The figure shows a series of snapshots of the 
disk's density in the rotational equator of the neutron star, which fills the central 
circle - the spatial scale is in units of computational cells, with 10 cells equaling one 
neutron star radius. The time sequence goes from left to right, with the lower row 
following the upper row. The simulation begins with the small disk of the equilib- 
rium electrosphere. At later times the disk spreads and develops nonaxisymmetric 
rolls and fingers, characteristic of Kelvin-Ifelmholtz instabilities, to which diocotron 
instability is closely related. 
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Fig. 7. Location of proposed charge starvation gaps in the magnetopshere, 
employed in models f or pu lsed high energy emis sion. Left: slot g ap model o f 
lArons fc ScharlemannI ll 19791) : IXrons ( 1983bl: Mu sUmov fc Hardind (|2003l . \2004 ): 
figure from Muslimov fc lardind I2OO3I'). Right: c l assica l outer gap geometry of 
ICheng Ho fc RudermanI l| 19861) : IChiang fc Romanil l| 19941 ) :[ Romani fc Yadigoroglul 



l l995l ): lRomanil (|l996l ): figure from iHirotanil l|2006t ) 



These models use magnetic geometry borrowed from the vacuum obUque 
rotator, with plausible but ad hoc prescriptions for the shape of the sepa- 
ratrix and for the choice of field lines assumed to participate in the gap. 
Strictly steady flow (and therefore electrostatic accelerating electric fields in 
the corotating frame) are assumed, with a variety of higher order drift effects 
included on the high energy particles' orbits, in order to capture the rather 
highly refined beaming dynamics at large radius, including pulse structure 
controlled by caustic formation through delicate cancellations of aberration 
and light travel time with field line s weep back along a last closed flux sur- 
face of assumed form (fPvks fc R,udakH20mi ). Manipulation of a variety of free 
parameters (especially the thickness of the assumed accelerating layer) allows 
fits of the resulting radiation spectra and pulse profiles to observations with 
greater or lesser success. Since such geometric constructions are sensitive to 
the exact form of the geometry guessed, a fruitful exercise would be to place 
such models, or alternates such as the gapless model involving the return cur- 
rent itself outlined below, in the context of the force-free solutions, where the 
location and shape of the separatrix is not a free function, or is a function 
only of Ry/Rl- At the very least, such model construction would be a step 
toward probing the basic structure of the magnetosphere, a task made possi- 
ble since all the phenomenological gap models contain parallel voltage drops 
Z\^ll — J E^\ds small compared to therefore allowing the force-free theory 
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to be a good zeroth order platform for parallel accelerator and radiation model 
construction. 



2.6 Current Flow Profile and Gap Electrodynamics 

The force free electrodynamic solutions also exhibit an important result 
which affects all the gap models construct ed for the last 30 years, s tartin g 
with the vacuum surface and o uter gaps of iRuderman fc Sutherland 



Cheng. Ho fc: RudermanI (|l986f ) and their many successors, as well as the 



space charge limited flow be a m models, both with and without slot gaps, of 
Arons fc ScharlemannI (|l979t) : lAroiii (|l981al Il983bl ) and the many successors 



of this modeling idea. 

All these schemes embody the idea that appears as a result of the 
magnetosphere's attempt to restore charge neutrality in the co-rotating frame 
by accelerating non-neutral beams of particles with density comparable to 
the Goldreich- Julian density rjji/q, with that adjustment to perfect charge 
neutrality in the corotating frame being incomplete - an due to charge 
starvation. In such configuations, the field aligned current J|| = J ■ B/B 
adjusts to a value controlled by the local electrostatics of the hypothesized 
region of charge starvation, not to a value determined by the energetically 
dominant magnetospheric dynamics. 

For currents emerging from the star's atmosphere, the accelerating E'n ap- 
pears from the gravitational depletion of density below the Goldreich- Julian 
value, due to the low temperature and low radiation pressure which prevents 
the filling of the magnetosphere and formation of a wind with a charge neutral 
plasma pushed up by pressure and centrifugal forces from the stellar surface. 
The current supplies the charge needed to support a polarization electric field 
which (almost) cancels the vacuum field. The residual ("starvation") electric 
field still supports enough of a voltage drop to allow the accelerated parti- 
cles to emit gamma rays that can convert to electrons and positrons. The 
resulting gap structure, elaborated assuming strictly steady fiow in the co- 
rotating frame, thus enforces a current density almost uniform with distance 
from the magnetic axis, with value close to the canonical value crjR. Since 
this current fills (almost) the whole polar cap, the total current from a polar 
cap is / « C77fl7rr^^p « c{Q B cap cos i / 2t: c)t: / Ry — c<P, which suggested 
that such a gap might be an element of the magnetospheric circuit, although 
with the peculiar property that the charge density (and therefore the current 
density) is an eigenfunction of the local electrostatics. The return current is 
not included in these local acceleration models, being explicitly or implicitly 
assigned to the current sheet. 

Such local determination of J|| is not what one expects on energetic 
grounds, since the current density refiects the induced magnetic field, through 
which all the spindown energy loss flows. That energy fiow is much larger 
than the energy fiowing through the proposed electrostatic accelerator, which 
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is thought of as a smaU perturbation of the force-free structure. In the ab- 
sence of further information from full magnetospheric solutions, or specific 
features coming from phenomenological models of radio or high frequency 
pulse obsevations which characterize the current flow in more detail, the hope 
expresssed by the approximate correspondence between the total gap current 
and the current of the magnetsopheric circuit has stood unchallengecQ. 

For outer gaps, whose morphology appears in Figure[7l the current density 
is also established by the effect of pair plasma shorting out the starvation , 
since within the gap counterstreaming electron and positron beams coexist, 
with acceleration ceasing at the end points where the pair density rises to be 
approximately equal to the Goldreich-Julian density - or not at all, along field 
lines closest to the (assumed) conducting boundary formed by the last closed 
field line, where pair creation is weak, for reasons traceable to the assumed 
geometry. Figure [5] shows the place outer gaps might have in a hypothetical 
picture of the global circuit. The starvation electric field can be sustained 
only if the beams have density not exceeding crjji/e, therefore the current 
density necessarily approximates the Goldreich-Julian current density. Since 
the gap must be thin in the poloidal direction across B (otherwise the photon 
emission from the gap would not make a narrow pulse, the express purpose 
for which the gap model was designed), the total gap current has to be small 
compared to the magnetospheric current - outer gaps cannot close the whole 
magnetopsheric circuit, if they are to have any pretensions as a successful 
model for pulsed gamma rays. Their biggest success, which they share with 
the slot gap model, is the assignment of the radiating geometry for gamma 
ray pulsars to a thin sheet which closely follows the last closed flux surface of 
the oblique rotator. 

Outer gaps, as regions of field aligned acceleration relying upon starva- 
tion electric fields, can occur only on field lines not supplied with a dense 
plasma either from the polar caps, as comes from the polar and slot gap space 
charge limited acceleration region, or from the recirculation of polar outflow 
by reconnection flows in the vicinity of the Y-line. Also, the outer gap, if it 
exists, sends almost all o f its pair plasm a back to the stellar surface, rather 
than supplying the wind ( Hirotani 2006[ ). Thus the outer gap has a hard time 



being a major supplier of the known large {N± ^ Nn) particle fluxes known 
to be injected into the young Pulsar Wind Nebulae, as discussed further in 
^ Outer gaps also run the risk of supplying too much energy in precipitat- 



® That one might be able to use observations to probe the current flow structure has 
be en an almost untouch e d subj ect. One of the few counterexamples was provided 
bv lHibschman fc Aronsl (|2001a!). who pointed out that thin return current layers 
might create observable signatures in the radio polarization data. Their predicted 
signature of the return current layer may have been seen in PSR J1022-I-1001 
(|Ramachandran fc Kramerll2003l ). 
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Fig. 8. Current flow patterns hypothesized for modern outer gap models (jHirotanil 
[2OO6. ) , which include possible outward ion flow on field lines near the separtrix (in 
i < 90° geometry). Left: global pattern assuming a closed circuit, with currents 
closing in the wind. Right: hypothetical current flow from the star through the 
nuU surface. No attempt is made to account for the return current in the current 
sheet in these models, nor has there been any correspondence made between these 
hypothetical current flows and the currents in the force free solutions. 



2.7 Gap Subversion: Non-Uniform Current Profiles 

All gap models (vacuum polar cap gaps, space charge limited flow polar cap 
and slot gaps, outer gaps) function as suppliers of plasma which come as close 
as possible to restoring the charge density to the Goldreich- Julian density r]R. 
They rely on starvation electric fields, since the gaps' charge densities fall 
below rjji (by a little or a lot, depending on the model). Since the voltage 
drops developed in the gaps are highly relativistic, such plasma takes the 
form of relativistic particle beams, with the resulting current density parallel 
to B being J|| = c?7fl(l — h), where /i <C 1 for polar cap and slot gaps (the 
stellar surface supplies a beam of charge density almost equal to 77^^, thus 
trapped particle backflow formed at the PFF is small), while in traditional 
outer gaps, /i w — 1, (the pair formation front at each end has to trap plasma 
with the full Goldreich- Julian density, forming two counterstreaming beams of 
approximately equal density). This current is almost constant, as a function 
of distance across the magnetic flux surfaces (formally, J\\{iIj) ~ constant). 
All the models assume strictly steady current and plasma flow in the co- 
rotating frame (electrostatic approximation), even though the models' authors 
frequently indulge in discussion of time dependence that they think should be 
part of their proposals. Steady current flow has been found in evolutionary 
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force-free and relativistic MHD models I Spitkovskv 2006t [Koniissarovll2006l) 
these show no signs of variabihty in the corotating frame (they also cannot 
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Fig. 9. Polar current flow patterns hypothesized and found from force free models. 
Top Left; Polar current flow with all return current in the current sheet, and illus- 
trating balancing the charge loss with an ion current extracte d from the surface, in 
response to electrons and plasma precipitating from the Y-line (|Aronj|l981bl . ll983al ) 
Top Right : Magn etic field lines of the aligned rotator when Ry/Rl = 1-0, from 
[McKinneyl (|2006l ). showing the asymptotically monopolar character of the poloidal 
field structure. Bottom; Current flow in a force fre e solut ion of the aligned rotator, 
for various values of 2:0 ~ Ry/Rl, from iTimokhirJ l|2006l ). 9 is the magnetic colati- 
tude of a field line's footpoint, and 9pc is the magnetic colatitude of the polar cap's 
edge, defined as the intersection of the separatrix between closed and open field lines 
with the stellar surface. 
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capture reconnection physics, which probably requires higher resolution simu- 
lations than have been employed so far, and probably also incorporating more 
physical models of dissipation and inertia. But, even though global and local 
theories both have stationary lighthouse behavior built in, the consequences 
of gap electrostatics for the current flow distribution are entirely inconsistent 
with the results of force free modeling, because of a serious mismatch between 
the current profiles found in the gap and in the global models. 

Figure [9] shows the poloidal current density, as a function of magnetic flux, 
exhibiting the fact that a fraction (~ 20%) of the return current flows on open 
field lines just within the boundary of the closed zone if Ry = Rl] the rest 
of the return current lies within the unresolved current sheet separating the 
closed from the open field lines. Thus, as far as current flow is concerned, 30 
years of research on the force free rotator can be summarized by saying that 
to within 20% accuracy, the current flow distribution of the force free dipole 
is that of the monopole, whose total current appears in ([5]) - as far as the open 
field lines are concerned, the dipolar magnetosphere is the monopole mapped 
onto a polar cap, in each hemisphere. The current density as a function of 
distance from the magnetic axis then is (ignoring the small piece of the return 
current required on open field lines) 



dl . ( V' \ ■ / ^ 



,2 



dtp V VcapJ V ^pcJ 

where jcj = criii{'uj^ = 0) = QBcap cosi/27r and tu* is the cylindrical distance 
of a field line's footpoint from the dipole axis. If the polar current is a charge 
separated, steadil}!^ flowing beam extracted from the star's atmosphere by 
E\\ with charge density 77 = J\\/cl3, with (3 k, 1 except in a thin region at 
the surface of thickness 10 or so times the atmospheric scale height, then the 
difference of this beam charge density from the Goldreich- Julian density is of 
order "q^ itself over most of the polar flux tube: 

■h_ ^ jcj L ijj \ jcj ^ ^ be .7 1 ^ \jGJ 



Vr = 1 ; — = H 5 > : ^ ^cap- (10) 

C C \ IpcapJ C C W^ap C 

Expression (fTU|) means that the parallel electric field is not almost shorted 
out by the space charge density of the particle beam from the surface, thus 
returning the space charge limited fiow to an environment with an electric 
field akin to what was envisaged for vacuum gaps (if the plasma forming 
the return current on the polar flux tube's boundary behaves as a perfect 
conductor, as has been assumed in all models to date), but now with the 
requirement (since the stars have dense, thermal X-ray emitting atmospheres) 
that E\\{r = i?,) = 0. 

The resulting hu ge acceleration inevitably leads to massive pair creation. 



in the manner of the lSturrockl (|l97ir ) picture of a high energy beam coexisting 



^ Steady on times long compared to the polar cap transit timescale zucap/c ~ Wfj, 
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with massive pair creation and acceleration with almost all of being dropped 
by the electric field within a height about equal to the polar cap width. Such a 
situation has both interesting possibilities and large problems, both theoretical 
and observational. 

Pairs forming in an approximately vacuum electric field will short out at 
a height such that the voltage drop is sufficient to allow the accelerating 
particles of the beam from the surface to emit magnetically convertible gamma 
rays. Typically, Zi^n ^ fO^^ V ^ ^. Because the difference charge density 
pI7|) is a large fraction of —rjji over most of the open flux tube, the pairs 
have to supply most of the charge density needed to shut down at the 
pair formation front (PFF). Because the pair formation front is now formed 
at lower altitude than is the case in space charge limited flow models with the 
current density determined by the local electrostatics, pair cascades may be 
better able to supply the total plasma flux inferred from pulsar wind nebulae, 
as discussed in ^21 However, the pairs generally are born with energies small 
compared to eZ\<?. Therefore, particles with charges having the same sign as 
the particles of the beam extracted from the stellar atmosphere are added 
to the outbound beam, while particles with the opposite sign of charge are 
trapped electrically and go backwards - backflow from the PFF enhances the 
current. If the potential is monotonic, the particle backflow collides with the 
stellar surface with number flux « {c\r]ii\/e){'Uij/wcap)^ and energy /particle 
eA<P. The particles in the backflow lose their energy after penetrating several 
hundred gm/cm^ into the crust, heating the atmosphere from below. Then 
each pole would have thermal X-ray emission with luminosity/pole Lxpoie ~ 
0.5Er{A$\i/$). With A<P\\ ^ 10^^ V (the voltage defining the theoretical 
death line in Figure [T]), these heated polar caps would emit substantially 
more thermal X-rays than are observed in many RPPs. Similar backflow in 
"spark" models runs into the same difhculty. 

The traditional space charge limited polar caps, in which the current car- 
rying beam extracted already has density close to |?7jj|/e, greatly reduce this 
emission - in modern models, in which dragging of inertial frames controls the 
difference between the beam's density and |?7i?,/e|, the reduction is by a fac- 
tor ~ OAGM/ R^c^ « 0.06. This reduction is enough to give polar cap X-ray 
emission either in accord with observations of some pulsars, or small enough 
to be hidden by magnetospheric non-thermal emission. But the price paid is 
an electric current (formed by response to backflow from the PFF as well as 
by emission from the atmosphere) over the whole polar cap (stationary or 
non-stationary) which is large compared to (O exce pt near the magnetic axis. 
Furthermore, just as in the early spark gap models of Ruderman fc Sutherland! 



( 1975f ). the difference charge density in (fTO|) implies a nonzero E\\ incident on 



the star's surface - really, the top layers of the atmosphere 

A number of ideas have been advanced to resolve this conundrum. 

• The PFF has a different structure than has been found in studies starting 
with Arons Sz Scharlemann l|l979l) all the way through recent work on 
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the full slot gap ( Muslimov &: Hardind 2003 , 2004 ) . If somehow electrons 



could be made to precipitate from pair plasma above t he PFF and ente r 
the current flow region, the current would be reduced (jTimokhin [looi). 



However, this is quite hard to achieve, the electric field below the PFF 
acts to expel such particles - they could not enter unless they already 
had energy ~ eA<^\\ . For curvature gamma ray emission generated pair 
cascades, the pairs are born with energies far below eZ\^|| , which make 
formation of a PFF that reduces J\\ impossible. 

The PFF and the acceleration region below it is fully time dependent, 
with the current flickering on the p olar cap t ransit time ^ vJcav/c 
10 Mse c, in the manner of the Sturrockl (|l97li ) and lRuderman fc Sutherland 



(|l975l ) diode instability and spark scenarios, respectively. Given stars with 
atmospheres which make charges freely available to the magnetospheres 
above, the most significant gap question is, how does the star adjust to 
provide the magnetospherically rerquired current density, while at the 
same time adjusting the charge density to reduce E\\ below its vacuum 
value in a manner that does not do violence to thermal X-ray emission 
(for example). 

All the models of relativistic field aligned acceleration considered to date 
have assumed unidirectional flow, or at most counterstreaming beams with 
no trapped particles in the acceleration zone - that is, monotonic acceler- 
ating potentials and electric fields. The possible presence of trapped parti- 
cles, implicit in the non- monotonic, non-relativistic acceleration model o f 
a singly charged fluid of IL. Mestel and M. Prvce. in Mestel et al. ( 1985[ ). 



can break the straight jacket of current being proportional to charge den- 
sity, even if the plasma is completely non-neutral. Furthermore, the outer 
magnetosphere can modify the local polar cap (and outer gap, if it exists) 
electric field through return current control of the currents and charge 
densities in the polar flux tube boundary layer separating the open from 
the closed field lines, which opens the possibility of relaxing the assump- 
tion of perfect conductivity on the boundaries of the acceleration region 
used in all models to date. All models to date have assumed the open field 
lines are bounded by a surface whose behavior mimics a perfect conductor, 
which makes a very specific statement about the surface charge density 
along the separatrix. As with the total current and current density, there 
is no reason for the surface charge to adopt this locally determined value 
- it depends not only on the charges contained within the open flux tube, 
but also on the dynamics controlling the formation of the return current, 
both at the stellar surface and at the magnetopsheric Y line. I discuss 
these speculative points further below. 

The magnetosphere is a high inductance system - the magnetic fields 
induced by particle and displacement current flow can change on times 
no shorter than the Alfven transit time between surface and Y line and 
back, Ta ~ 2RyIva ~ {'2'/ ^){Ry / Rl), very much greater than the po- 
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lar cap transit time, or the plasma period based on the Goldreich-Juhan 
density uOpQj — A'n:e^\rjpi\/m±e — ^2Qflc± — fl\j2ftc±l fl ~> fi [here 
flc± = eB/m±c is the nonrelativistic cyclotron frequency ~ 10^^(i?/10^'^G) 
of an electron or positron]. Thus one possibility is that the atmosphere 
supplies the current demanded through pair creation discharges, which launch 
charge bunches accompanied by pulses of much denser pair plasma, a popular 
idea frequently mentioned in the cartoon approximation but rarely studied. 

Homogeneous and spatially ID mod el s of time dep ende nt pair creation 
have been studied bv lAl'ber et'aR (Il975l ). iFawlevI and iLevinson et al. 



(|2005^ - the first considers only time dependent oscillations in a uniform 
medium (a "OD" model), the second, using a ID Particle- in-Cell method, con- 
siders the propagating transitions between vacuum and an ^ region 
creat ed by a burst of pair creation (motivated by the Ruderman fc Sutherland! 
19751 spark scenario), while t he third cons i ders n onlinear uniform oscillations 



similar to those studied by lAl'ber et al\ (Il975f) and also spat ially inhomo- 
geneous nonlinear limit cycle oscillations. Levinson et al. ( 2005( ) incorporated 
the existence of the magnetospheric current as a fixed constant in the model, as 
is appropriate since the oscillations occur on the time scale i^~qj, much shorter 
than all possible time scales of magnetospheric variability. They point to the 
interesting possibility that the charge and current oscillations might become 
chaotic, a topic of substantial interest to possible radio emission mechanisms, 
but present no specific calculations that exhibit such behavior. 

None of these explicitly time dependent models showed approach to the 
steady flow in the corqtatin g fra me assumed at the start in the models of 
Arons fc ScharlemannI (|l979l ) and iMuslimov fc Hard"in3 (|2004 ). in which the 



current density is fixed by the local electrostatics. Significantly, none of these 
local, OD or ID time dependent models included the effects of pair outflow 
from the system, nor included the Poynting fluxes (in effect, the collective 
radiation losses) from the regions of pair oscillation. Considered as an insta- 
bility, the oscillations studied are probably sensitive to the loss of plasma 
and Poynting flux from the system, since the pair creation (in ID) has the 
character of a spatial amplifier. Thus the question of whether pair oscillations 
approach a steady state with the local value of the magnetospheric current 
density emerging as a steady flow or as an average over the oscillations re- 
mains open, requiring calculations which are able to give a more complete 
account of the coupling to the stellar surface and the losses from the region 
of oscillatory dynamics. Giving an account of the effects of plasma flow onto 
the star and the consequent heating and X-ray emission is of substantial im- 
portance, since the Levinson et al. model suggests a local current density crjn 
accelerated through ~ 10^'^ V colliding with the surface, which leads to ob- 
servable thermal X-ray emission which may, or may not, be in excess of what 
is observed from many sta rs. 

On the opposite side, lArond (|l983bl) and iMuslimov fc Hardind (Il997l) 
point to how their steady flow, spatially inhomogenous models ("gap-PFF") 
might become unstable, due to inhomogeneity of the pair creation, in older 
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stars where the pair creation gain lengths are large. Such instability would 
be of the spatial, traveling wave amplifier variety, a possibility also of inter- 
est to the outer gap models, whose local time dependence has also begun to 
be studied. iHsu. Hirotani fc Chand (l2006l) have opened the first door to time 



dependent (in the co-rotating frame) outer gap models, showing that their 
models, which do incorporate particle outflow from the accelerator region - 
mostly toward the star - converge to a steady state flow. Strictly speaking, 
however, these models are inconsistent, allowing for full time dependence of 
the current and plasma densities but treating the electric field as electrostatic, 
which is quantitatively incorrect in a relativistic, multi-dimensional (2D, in 
their model) system. The approach to a steady state is attributed to the 
screening of if too many pairs materialize, while underscreening results 
in an increrase of E\\. This mechanism is the same as underlies the nonlin- 
ear limit cycle oscillations appearing in the polar cap pair oscillation model 
of Levinson et at. Thus the approach to a steady state in the time depen- 
dent outer gap model more likely owes its origin to the spatial loss of plasma 
from the acceleration region, an effect broadly akin to transit time damping 
of electromagnetic oscillations in a plasma. 

Given the large, almost vacuum conditions above much of the polar cap 
implied by pH]) . coupling to the stellar atmosphere almost certainly requires 
consideration of trapped particles' contribution to the charge density, in either 
steady flow or time dependent current flow conditions. 

A full theoretical resolution of pair creation driven oscillations in any "gap" 
geometry awaits more definitive study, along with any sort of serious attempt 
to relate such oscillations to observable phenomena - current and torque fluc- 
tuations, radio microstructure, variability in thermal X-rays created by surface 
bombardment, etc. All such modeling needs to be set into the context of the 
global force-free models - so long as the potential drop in a local accelerator 
is small compared to the magnetospheric potential fleld aligned ac- 
celerators of any sort (employing starvation E^^ in all the schemes available in 
the literature) appear as small departures from force- free conditions, allowing 
the use of force-free models as the basic zeroth order description of the current 
flow and magnetic geometry. 

2.8 Gamma Ray Tests of Existing Gap Models 

The modern force-free magnetsophere models o pen the possibility of usin 



the upcomin g gamma ray observatory GLAST (jGehrels fc Michelsonlll99' 



Carson 20061 ) to test and improve our understanding of pulsars' magneto- 



spheres, along the lines suggested above or in other directions. Gamma rays 
afford the possibility of probing the magnetosphere using well understood ra- 
diation processes, leaving the modeling and the synergy between models and 
observations living in the domain of the geometry and the acceleration physics 
- "gapology", in the existing theoretical frameworks. In particular, the advent 
of the force free-models should allow the outer gap and slot gap modelers 
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to significantly reduce the geometric uncertainties in their constructions of 
the beaming profiles and energy dependent light curves, thus allowing much 
more stringent empirical tests of the hypothesized geometric scenarios - if the 
necessary extensions of the models to 3D, and perhaps to time dependence, 
are incorporated. Even more important, the improved sensitivity of GLAST 
over past gamma ray telescopes will allow, for the first time, a direct test of 
whether polar caps and polar cap pair creation occur in a significant pop- 
ulation of pulsars. Given that no proposed outer gap (or slot gap, for that 
matter) model makes a significant contribution to pair creation and gamma 
ray emission for periods much in excess of ~ 200 msec < 10^^ Volts), the 
much heralded association of pair creation with pulsar photon emission and, 
more significantly, with relativistic wind formation must come from activity 
in the polar cap region just above the surface. Indirect evidence for such pair 
creation comes from the simple 'P = 10^^ V radio pulsar death line shown in 
Figure [H which corresponds roughly to wher e polar cap/slot gap acceleration 
models predict pair creation to cea se (jSturrock 1971: Ruderman fc Sutherland 
1975t lArons fc Scharlemann 1979f l. Previous gamma ray telescopes lacked the 
sensitivity to probe the predicted gamma ray emission, which, in the models, 
is absorbed at energies above 1 GeV (in simple, star centered static dipole 
geometry) through gamma ray conversion to pairs. 




Fig. 10. Sensitivity of the EGRET and GLAST telescopes to low altitude (below 
and just above the PFF) gamma ray emission at energies e > 100 MeV, in a simple 
dipole model for the low altitude geometry (which affects the magnetic curvature 
and therefore the maximum photon energy that can escape) and unidirectional space 
charge hmited flow with current density ~ cr/R, not the force- free current given 
(approximately) by Q. The different symbols re fer to the different major contribu- 
tors to gamma ray emission and absorption - see iHibschmanl (|2002l ') for the details. 



Pulsars; Progress, Problems and Prospects 27 



Figure [TOl taken from iHibschmanl (j2002f ) , gives a simple version of this 
opportunity, in good accord with more recent evaluations of pulsar gamma 
ray emission from the inner magnetosphere. Testing the existing polar gap 
acceleration and gamma ray emission models, or better, improved models 
that take proper account of the magnetospheric current system, can be done 
best by studying pulsars that show core component radio emission, since one 
looks down the "barrel of th e gun" into the low altitude magnetic field, where 
the core emission arises (see iKramer et al. I ll99al for the evidence that the core 
emission arises from a few kilometers above the surface in a substantially 
dipolar B field). 

I can safely predict that GLAST observations of pulsars deeper in the 
P — P plane than could be detected by previous gamma ray telescopes, which 
should begin to become available in 2008-09, will stimulate a small host of 
improved gap models which take advantage of the geometric and current flow 
constraints coming from the force-free models. They might stimulate investi- 
gations into origins of i^y based on physics differing from the starvation models 
that have been explored for the last 30 years - for example, invoking the -Ey 
accompanying the kinetic Alfven waves that couple the outer magnetosphere 
to the star, a conceivable acceleration mechanism that might augment or even 
replace outer gap and slot gap accelerators, especially if plasma precipitating 
from reconnection flows at the Y-line floods the regions envisaged for these 
gaps and poison their starvation electric fields. 

To amplify this issue, which is a prospect for future research, consider the 
hypothetical global reconnection flow illustrated in Figure [HJ shown here for 
an aligned rotator but just as applicable to the Y-line of the oblique rotator 
illustrated in Figure [3l 
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As appears in the iBucciantini et model, pairs supplied from 



the polar gap supplying the wind should allow reconnection of the current 
sheet to occur all the time. It has not yet been shown that the reconnection 
propagates back into the magnetosphere in the manner envisaged in Figure 
[TTl In the somewhat analogous problem o f reconnec t ion at the Y line in the 
rapidly rotating Jovian mangnetosphere, Yin et al. ( 2000f ) showed, using a 



PIC simulation of the electron-ion plasma at the Y-line, that reconnection 
qualitatively akin to that illustrated in Figure H] occurs, with finite but not 
large reaction along the separatrix interior to the sporadically forming X- 
line. That back reaction includes generating precipitating Jy and plasma with 
density well in excess |J|||/ec on and around the separatrix. The field aligned 
currents are part of the kinetic Alfven waves that couple the time variable Y- 
line to the inner magnetosphere and the star, thus generating a time variable 
torque. The space charge in these boundary layer flows can alter, in a major 
way, the electric field within the polar cap accelerator (as well as poison both 
outer and slot gaps), while the parallel electric fields in the kinetic Alfven 
waves offer a new mechanism for field aligned acceleration in the boundary 
layer geometry already known, from the outer and slot gap models, to be 
extremely useful in understanding gamma ray pulsars' beaming morphology. 
If these waves have a chaotic time series, the chaos in the resulting polar Jy 
offers a good candidate for understanding the random arrival times of radio 
subpulses within a pulse window; if the reconnection induced waves have a 
limit cycle time series, the phenomenon of subpulse drifting can be reproduced, 
assuming the radio emission intensity and beaming is a direct product of the 
field aligned current density. The electric fields in these Alfven waves have a 
central role in the formation and extraction of the return current required to 
maintain the average charge balance of the star. Finally, the shear between 
the plasma flow in the polar flux tube and the neighboring boundary layer 
offers a promising candidate for the collective radiation mechanism in "conal" 
radio emission. 

An important constraint which must be met by any model of outer mag- 
netospheric variability is that magnetic fluctuations must not broaden the 
beaming of photons emitted with momenta parallel to the instantaneous mag- 
netic field beyond the characteristic pulse widths observed in the gamma-ray 
and associated optical and X-ray em ission (a.nd rad io emission, in the Crab) , 
thought to come from r > 0.5Rl (jRomani Il996^ - although this number, 
derived from a geometry based on a vacuum dipole with a phenomenological 
prescription for the location of the last closed field lines, will change when 
the magnetic field of the oblique force free rotator is put to use as the geo- 
metric platform for the beaming. If reconnection (or any other mechanism) 
causes the magnetic field at the light cylinder to fluctuate by an amount SB at 
r Rl, then the Alfven waves traveling along the boundary layer have, from 
conservation of energy flux, amplitude 6B = B{Ri,){8B / B)r=R^{r / RlY/"^ - 
for convenience, I have here assumed Ry = Rl- These are shear waves, with 
6B _L B, thus causing the local magnetic direction to vary about the mean 
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by an angle 69 w SB/B. The observed sharpness of the gamma-ray pulsars' 

Hght curves then suggests 59 < 0.2{B /SB)^^Jj^^, while order of magnitude 
application of this idea to the observed torq ue fluctuations in the Crab pulsar 
and others suggests [B / 5B)r=RL < 2 ( Aronsl[l981b ). The correspondence of 
the limits on the emitting radius from geometric fluctuation pulse broadening, 
from measurements of the torques and from modeling the beaming geometry 
provides an example of how gamma ray observations of pulse profiles and ra- 
dio observations of torque variability over a substantial range of the P — P 
diagram can be used to seriously constrain both the electrodynamics and ac- 
celeration physics of these magnetospheres. 

I have left out all discussion of the hoary problem of pulsar radio emission 
and transfer - that would require an additional paper - other than the few 
comments above concerning radiation beaming and single pulse fluctuations, 
which appear to me to provide probes into the magnetospheric dynamics. In 
connection with the dynamical importance of the boundary layer bedween 
closed and open field lines - the location of the return current - it is perhaps 
worth emphasizing that this region is likely to be the dynamical realization of 
the site of "conal" emission, with velocity shear between the boundary layer 
plasma and the plasma filling the open flux tube (and that filling the closed 
zone) as a prime candidate for the free energy driving the collective radio 
emission process(es). 

Perhaps a few bold souls will explore these issues more quantitatively 
before the GLAST observations become available, thus offering up their pre- 
dictions to the sharp knives of experimental tests. There is an urgent need 
for physical models of the boundary layer between the open and closed re- 
gions, either with or without gaps, which account for the coupling with the 
stellar surface as well as the transition from the magnetosphere to the wind. 
This is a collection of non-trivial problems - predictions of future progress are 
uncertain. 



3 Follow the Mass 

While there has been lots of attention to pair creation and particle acceler- 
ation within pulsars' magnetopsheres from the community interested in ob- 
serving and modeling these stars' lightcurves and SEDs, the most obvious 
evidence for particle particle acceleration and pair creation comes from ob- 
servations and models of Pulsar Wind Nebulae (P WNe) . T hese h a ve rec ently 
been well reviewed from the observational standpoint by Island (|2005h and 



Gaensler fc Slane ( 2006f l. The young PWNe and their pulsars - those still not 



crushed by the reverse shock - provide calorimetric information on both the 
energy and mass loss budgets of the underlying pulsars. Indeed, since the ear- 
liest days, the energy budget has been used to constrain the moment of inertia 
of the neutron stars, thus the equation of state of nuclear matter. The mass 
loss budget provides a powerful constraint on all models of plasma behavior 
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within the magnetosphere, whether designed to explain specific observations 
or constructed to investigate basic theoretical issues, 



3.1 Observations and Consequences 



There has been major observational progress on these systems, coming most 
of all from high resolution optical and X-ray imaging, as shown in Figure 1121 
and from related studies of temporal variability. 

Observations and models of the PWNe tell us about the particle loss rates 
7V-I-, Nheam from the neutron stars. In the case of the Crab, with its strong 
magnetic field and rapid synchrotron cooling of the particles radiating photons 
at energies above infrared frequencies, the now well resolved optical. X-ray 
and the unresolved gamma ray sources require a particle input of around 
-[^q38.5 g-i^ about 1 0^ times the basic electric current flow c<P/e ~ lO'^'*'^ s~^ 
for this pulsar (e.g. 'Spitko vskv fc Arond 2004 ). Similar conclusions have been 
reached for oth er PWNe as has been done, for example, in G320 around 



PSR B1509-58 (jGaensler et aU 120021: iDeLanev et al\\20m . even though in 



this case the radiative l osses are n ot as rapid and therefore inferring N is not 
as straightforward. See Ide Jagei ( 2007 ) for pair injection rate inferences for 
several other pulsars/PWNe. These inferred rates come from examining the 
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Fig. 12. Left: composite of Hubble Space Telescope and Chandra images of the 
inner 1' of the Crab Nebula, showing the torus-plume structure around (torus) and 
along (plumes) the pulsar's rotation axis. Right: Chandra image of the supernova 
remnant MSH 15-52 (G320.4-I-1.2), showing its torus-plume structure. 



brightness of the X-ray and (when they can be seen) optical nebulae, aug- 
mented by hard X-ray and gamma ray (GeV to TeV) observations of the 
young nebulae, when the hard photon telescopes can detect anything - for a 
recent summary of the rapidly developing TeV observations of TeV PWNe, 
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see I Gallant) (j2007f) and iFunld (j2007f ). The TeV observations are sensitive to 
particles of energy comparable to those that give rise to synchrotron X-ray 
emission. 10 MeV to 10 GeV observations constrain the highest energy par- 
ticles, which produce synchrotron X- and 7-ray emission. A survey of PWN 
emission in this energy range awaits GLAST. 

While the pair production rates found from polar cap models based on 
starvation electric fields do seem a dequate for the understandi ng of the high 
energy photon emission in PWNe ( Hibschman fc Arond 2001c - these results 
remain the only attempt to survey pair creation across all observed P — P), 
other models have been developed specifically for the purpose of explaining 
gamma ray pulsar SEDs and light curves), there are substantial indications 
that somethin g of qualitative significance is missing. It has long been known 
{e.g. IShklovsk v 1968) that the total radiating particle content of the Crab 
Nebula (mostly in the form of radio emitting electrons or pairs, which lose en- 
ergy only because of adiabatic expansion) requires an injection rate averaged 



over the 1000 year history of the system on the order of 10 
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10 



,41 



elec- 



trons plus positr ons, in order to u nderstand the total radio emission from the 
Nebula. Recentlv lde Jagej ( 2007 ) has revisited this same question in the fight 
of the TeV observations of VelaX, G320 and the newly discovered nebula of 
PSR B1823-13, again finding pair injection rates greatly in excess of the rates 
found for particle outflows from either polar cap/slot gap or outer gap modelarl 
constructed using starvation electri c fields shorted ou t by the pair creatioio- 
Only the early polar cap model by iTademarul (|l973h . in which the effect of 
pairs' ability to limit the voltage drop was completely neglected, comes even 
close to yielding the observed time average injection rates. Since incorporating 
the pairs' polarizability destroys Tademaru's model, his empirical success has 
been ignored. 

The starvation electric field models also have difficulty in coming up with 
enough pair plasma to meet the desires of most (not all) models of rad i o emis - 
sion over the whole PP diagram. The results of iHibschman fc Arond (2001b) 



However, his conclusion requires extrapolation of the particle spectra inferred 
from the TeV emission to radio emitting energies, a big jump. 
Some outer gap models applied to the Crab pulsar do find total particle pro- 
duction rates in the range 10^^ — 10^" s"'^, but these refer to particles flowing 
in toward the st ar, where they c ollide with and are absorbed by the surface and 
(over)heat it. See lHirotanil (|2006l ) for the most recent version of this kind of model. 
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appearing in Figure [T3] show clearly that for lower voltages {<P < 10 
where most pulsars lie, the pair multiplicities drop well below unity, 
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Fig. 13. Multiplicities (num ber of pairs per particle in the Goldreich- Julian dens ity) 
across the PP diagram, from lHibschman fc Aron"^ (|2001bl ). lHarding et al\ (|2002l 'l re- 
port similar results, using a more elaborate analysis. Both used similar versions of 
space charge limited beam acceleration in the polar cap region, and both included 
the contribution of synchrotron cascades to the total multiplicity. Crosses refer to 
objects where curvature emission provides the gamma rays that convert to pairs, 
circles to objects where non-resonant inverse Compton upscatter of thermal X-rays 
(both from the polar cap heated by backflow bombardment and from the whole 
surface of he cooling neutron star) provides the gamma rays, while asterisks show 
the more strongly magnetized stars where the cyclotron resonance in the scattering 
cross section makes a significant contribution to the gamma ray production rate. All 
these calculations used a star centered dipole for the magnetic geometry, and can 
be substantially affected by surface magneti c anoma lies, e.g. offsets of the dipole 
center from the stellar center, as in lAronsI (|l998bl ), or higher order multipoles, 
as in 'Melikidzc et al\ l|2006l ). The modern force- free model of the oblique rotator 
ISpitkovsky 2006) offers the possibility of investigating the pair multiplicity within 
a self-conistent geometric setting, either with the traditional starvation electric fields 
or with improvements that take into account the full magnetospheric current system 
and charge densities, as outlined in 



far below the level assumed in almost all models of the radio emission, and 
also required in models of propagation effects that have had some success in 
the interpr etation of rad io pol arization and beaming structure (e.g. McKinno3 
1997t iBamard fc Aro^lTsisl . or indeed needed to explain the death line in 



Figure [T] 
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It is possi ble that off s et of the dipole from the stellar center, a central 
aspect of the Ruderman ( 20061 ) model for magnetic field evolution, might 
substantially enhance a polar cap's pair yield, while still remaining consistent 
with the apparently dipolar morphology of the low altitude magnetic field. If 
the axis of the offset dipole is tipped with respect to a radius vector, gravita- 
tional bending of the gamma rays' orbits leads to much larger one photon pair 
creatio n, opac it y in t he magnetic field than is the case for the star centered 
dipole (jAronjll998bh . That opacity increase allows the more numerous low 
energy curvature gamma rays to contribute to pair creation, thus enhancing 
the particle fiux. This effect certainly does have a favorable impact on recon- 
ciling theoretical with observational death lines; whether it seriously enhances 
the pair yield in pulsars feeding plasma into young PWNe, which are far from 
death valley, remains an unexplored topic. 

The fact that pulsars (at least the young ones) must supply a plasma 
with particle outflow rate well above the Goldreich- Julian rate c^/e is un- 
deniable, based on the behavior of the PWNe - undeniable progress. The 
problems described here are quantitative, perhaps to be resolved by possi- 
ble larger pair output if the polar flow is time dependent ("flicker" flow), 
possibly to be resolved in a steady flow model if the charge density in the 
return current alters and increases the local electric field, perhaps by other 
effects not yet investigated. These possibilities have potential connections to 
time variability in the radio emission. Flickering of the polar cu rrent and pair 



creat i on might be connected to the radio microstr ucture {e.g. iJessner et al. 



iRuderman fc Sutherland ! 1197,4 lBenfordlll977h : modification of the elec 



trie field by boundary layer space charge controlled by outer magnetsopheric 
unsteady reconnection (see above) may be connected to subpulse variability. 
Improved sensitivity in infrared and shorter wavelength detection techniques 
that would allow probing for variability (in the corotating frame) of the higher 
frequency emission would be invaluable. Progress in this area in the next few 
years is to be expected. 



3.2 Pulsar Wind Nebula Models 
MHD Nebular Models 

Modeling the PWNe themselves has advanced greatly in the last decade. 
Drive n by the we alth of spatially and temporally resolved X-ray observa- 
tions ( Sland 2005 ) of the "torus-jet" structures shown in Figure [12] in the 
Crab and PSR B1509 nebulae and now known to be present in an increas- 
ing number of PWNe ( Ng fc Romanil 2004; Romani et al. 2005 ). modeling 
and si mulation has advanc ed fr om the elementary "sphe rical cow" mod- 
els of iRees fc GunnI (|l974l ) and iKennel fc Coronitil (|l984r ) to two dimen- 
sional, axisymmetric time dependent relativstic MHD simulations of the 
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flow structiire (IKomissarov fc Lvubarskvl 12003. 2004 iDel Zanna et a/.l 12004 : 
Bogovalov et aLll2005l) . 
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Fig. 14. Left: Outflow structure in MHD models with energy injection concentrated 
in the equator. A: Upstream relativistic wind. B: Subsonic equatorial outflow, down- 
stream of the equatorial termination shock. C: Fast downstream outflow emerging 
from the higher latitude oblique shock. D: Supersonic flow just downstream of the 
high latitude oblique shock. The backflow that foc uses downstream plasma onto 
the axis is not shown. F rom Del Zanna et al] l|2004l ) Right: Synthetic torus-plume 
image, from Zanna et aL (,200i ) 



Th ese simulations exploi t the suggestions of iBogovalov fc Khangoulian 
(|2002[ ) and lLvubars"k^ (|2002l ). that energy injected into these nebulae follows 
the cos^ A profile (A = latitude with respect to a star's rotational equator) of 
the toroidal fiel d energy density exhibited by the split monopole and oblique 
split monopole (|Bogovalovlll999l) models of the Poynting fluxes from the neu- 
tron star. They suggested that the consequence of such anisotropic energy 
injection into the surrounding nebulae would be greatly enhanced emission 
in a belt around the rotational equ ator - the " torus " appearing in the X- 
ray and optical imagery. In addition, Lvubarskvl 1 2002[ ) suggested the outflow 
from the torus, since it is injected into the nonrelativistically expanding cav- 
ity formed by the supernova, would be deflected into a subsonic backflow at 
higher latitudes, where magnetic hoop stress could act to focus plasma into a 
magnetically compressed, outflowing, subsonic plume along the pulsar's rota- 
tion axis, thus creating the appearance of a jet. 

The simulations amply confirm the implications of the initial toy models, 
with flow velocities in the equator and the plume {v ^ 0.5c) comparable to 
those inferred from motions of features in the published time series of motions 
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in the central regions of the Crab Nebula ([Hester et Snapshots taken 

from that time series are shown in Figure 1151 
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Fig. 15. The torus and plume i n the ce n ter of the Crab Nebula, as seen by HST 
and Chandra in 2001, from Hester et al\ |20o2). Left column; HST; right column; 
Chandra. Upper row; total structure, showing the knotty inner ring in the Chandra 
image morphologically identified with the wind's termination shock in the rotational 
equator. The elliptical shape reflects the 60° angle between the equatorial plane 
(which contains the torus) and the plane of the sky. Lower row; large scale structure 
subtracted by difi'erencing between a a pair of early and later images in the series, 
producing a snapshot of the moving wisps emerging from the Chandra ring and 
expanding at speed ~ 0.5c 



The (T Problem 

The MHD models do well at reproducing the torus-plume structure, as is 
shown in Figure [T4l if the wind upstream from the termination shock is weakly 
magnetized - the average of cr = B"^ /Airp'^c^ over the equatorial shock must be 
0.02 in order to generate the good looking simulated image, a value ^ 4x 
larger th an what had been p reviously inferred from the ID time stationary 



models of IKennel fc CoronitJ ([1984). But this value is still far below what the 
asymptotic cr expected in ideal MHD outflow of an unconfined wind, exhibited 
in (jlllP"! . and the asymptotic Lorentz factor and four velocity (in units of c) 



'^^ Relativistic MHD jets accelerating within confining boundaries - "walls" - 
may have quite differe nt behavior, as has been most recently exhibited by 
iKomissarov et all l|2007l ). Here, the effect of a confining boundary, perhaps rep- 
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2/3 ^ 1/3 /-i-ix 



2meffC^'yo 



leff ^ rriion + 2m±K±. (12) 



Here "f's = = Rl^, k± the pair multiplicity and 70 the bulk four 

velocity of the plasma emerging from the plasma source (pair creation at 
the polar cap, since outer gaps - if they exist as particle production zones - 
send most of their plasma back toward the neutron star) , itself an unknown 
function of magnetic latitude (possibly lower in the equator than at the poles, 
since pair creation should be weaker near the magnetospheric boundary layer). 
For the Crab pulsar, pair creation theory suggests the multiplicity of the 
plasma feeding the nebular op tical and X-ray source is k±ox ~ 10^-^, based 



on spectral imaging modeling (|Del Zanna et al. 2006) , a value consisten t with 
starvation gap modeling of the polar cap faibschman fc Arond 2001c ). who 



also find 70 ^ 100 for this star. If this piece of the mass loss budget corresponds 

to the total mass loss, Uo ~ 1 x 10'' and Wqo = 22 « 700. 

If one includes the whole mass flow, N± ~ lO""^'^ (jde Jageii 120071 ). 



which includes the particles required to feed the Nebular radio emission, then 
K± ~ 10® and (jQ is ^ 1Q'^^^ /jq - since the origin of the large mass flux is 
unknown, 70 is also unknown, although surely it is smaller than the value 
^ 100 found in existing gap models - then Uoc = 16(10/7o)^^'^. 

The ideal MHD values of a and 700 are for a wind with monopolar poloidal 
field and flow geometry at large r. Theory and simulation to date all support 
the poloidal field of the wind being monopolar well outside the light cylin- 
der - see, for example, iBucciantini et all (j2006t ). whose high a simulation of 
outflow from the aligned rotator extended to r « 900i?L, well outside the 
fast magnetosonic surface, with the field becoming closely monopolar with no 
polar focusing or hoop stress apparent. Then the asymptotic magnetization 
(Too and 4- velocity C700 (outside the current sheet, which is infinitely thin in 
ideal MHD) are predicted to be as in pTjl . 

The MHD models do answer the oft repeated question of just what is going 
on at higher latitudes - if the torus structure is the manifestation of the shock 
terminati on of the wind, w hy don't we see evidence for the shock at high lati- 
tude re.g.. lBlandfordll2002h ? The MHD model asserts that the polar regions at 



distances we can resolve are occupied by the backflow that forms the plume. 



resenting confining pressure in an outflow from a disk, forces poloidal field lines 
and stream lines to depart from monopolar form, which allows magnetic hoop 
stress to confine and accelerate a polar fiow. For a wind emerging from a star, 
essentially a point source, there is no analogue of confining walls to break the 
balance between magnetic hoop stress and electric repulsion that lies behind the 
slow acceleration. 
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The shock shown in Figure [TJ] curves down toward the star, reaching into radii 
too smah to observationally resolve on the polar axis. Furthermore, the shock 
is quite oblique at higher latitude, which weakens the efficacy of shock accel- 
eration. The MHD model and the curvature of the shock relies on the total 
energy flux being proportional to cos^ A. In the Crab, the higher latitude parts 
of the curved shock do manifest themselves as the bright knots, which appear 
in projectio n as if they are right next to the pulsar (jKomissarov fc Lvubarskv 



2003. 20041 ). Thus, qualitatively and semi-quantitatively, a satisfactory picture 



of PWNe plasma flow on nebular scales has appeared. 

Ideal MHD models may also do well at reproducing the wisp structures 
shown in Figure [TSl which are of interest for the diagnosis of the relativistic 
shock wave terminating the outflow. These are now known to be structures 
(probably waves, Scargic 1969) appearing to be emitted from the Chandra 
ring with a periodicity ~ 6 months, traveling out wi th a speed ^ 0.5c , pos- 



sibly with some deceleration with increasing radius ( Mori et al. 2006f) . The 



wisps occur on scales too small to be resolved by the published MHD simu- 
lations of the whole nebular flow. Of the various suggestions made over the 
years to interpret the wisps, the most promising MHD model for these is that 
they are due to MHD Kelvin- Helmholtz instabilities occurin g at the boundar y 
between a fast equatorial and a slower high latitude flow ( Begelma"nlll999() . 
The global flow models have such a shear layer, as flow emerging from the 
equa torial shock in the nebula r cavit y returns at higher latitudes toward the 
axis. iBucciantini fc del Zannal ( 2006h used high resolution MHD simulations 
of a shear layer in a box, repeated with periodic boundary conditions, to draw 
the conclusion that this hypothesis is quantitatively inadequate to reproduce 
the observed variability. However, recent (summer of 2007) high resolution 
MHD simulations of the whole nebular flow by Bucciantini and by Komis- 
sarov, both still in progress, suggest that either Kelvin-Helmholtz instability 
or, perhaps more hkely, secondary instability of the vortices formed in the 
shear layer near the base of the plume, do show many features similar to the 
observed moving wisp structures. 

If these models do exhibit as much ability to reproduce the observa- 
tions as has been found in the preliminary work, the multi-dimensional MHD 
model of PWN structure will have accumulated three major successes: 1) 
a model for thermal filament formation, through Rayleigh- Taylor instabil- 
ity of the boundary between the nonthermal bubb le, as first pointed out b y 
Hester et al. (1996) and developed extensively bv IBucciantini e^"a/] (12004): 



mildly relativistic plume (a.k.a. jet) fo rmation, as first suggested bv Lvubarskv! 
(l2002l) and mod e led nu meri cally by Komissarov fc Lvubarskv ( 2003. 2004h . 



Del Zanna et al. (|2004l ) and iBogovalov eTaLl (|2005h : and now the wisp vari 



ability near the termination shock. Such models probably will trun out to be 
successful in interpreting the more slowly expanding outer structures of the 
torus - the current round of high resolution simulations will soon show whether 
or not these featurs of the nebular "weather" can be captured numerically. 
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The MHD dynamics does have strong sensitivity to (Jwind at the termina- 
tion shock. The models are insensitive to the wind's 4- velocity (Lorentz factor) 
just upstream of the termination shock, and are insensitive to the composi- 
tion, other than that the particles must have small Larmor radii and that they 
be efficient radiators. The last requirement leaves electron- positron plasma 
as the only option, in the young systems with bright PWNe - the particle 
injection rates greatly exceed the Goldreich- Julian value. 



3.3 Beyond MHD 
Striped Winds 

However, the equator where the equatorial shock forms is a current sheet, a 
region notorious for breakdown of ideal MHD. Such breakdown has been as- 
sumed in the MHD models, which achieve their successful fits of nebular ap- 
pearance to observation only when there is a finite region around the equator 
where the magnetic field at the shock is small compared to what one would ex- 



pect in the ideal MHD flow with an infinitely thin, flat current sheet. iCoroniti 



suggested the apparent low value of a in the equator - infer red to be 



0.005 in the ID, spherically svmmetric iKennel fc Coroniti (1984) model 



is due to annilhilation of the equatorial magnetic field in the current sheet. 
Close to the star but outside the light cylinder the magnetic field takes the 
form of the striped magnetic structure, with oppositely directed fields from 
the opposite poles of the dipole wound into a frozen in wave, shown in Figure 

Coroniti's idea was that some form of current sheet dissipation causes 
the magnetic energy to annihilate in the inner wind, causing conversion of 
magnetic energy to flow energy, and reducing the structure to something ap- 
proximating that of an aligned rotator's outflow with a magnetic field in the 
asymptotic wind Rl ^ r <C Rshock (= 10^ Rl in the case of the Crab) much 
weaker than that what one expects from ideal MHD transport of the light 
cylinder field inferred from the star's spindown. The resulting (dissipative) 
MHD model has an equatorial current sheet built in, since the A > hemi- 
sphere has a toroidal magnetic field wound in the opposite direction to that 
found for A < 0. Figure [T71 shows a cartoon of the resulting magnetic "sand- 
wich" wind at large radius, along with magnetic field strength as a function 
of A considered in the MHD models of the nebulae beyond the shock. 

One almost model independent constraint on this idea is that an acceptable 
theory of stripe dissipation in the wind zone necessarily leads to the wind's four 
velocity in the dissipation region being small compared to the value Fwind 
IqG.5 inf erred from ID dynamica. 1 models of the nebular high energy photon 



spectra ([Kennel fc Coronitilll984l ). The reason is simple. 

The magnetic field in the stripes, which have proper wavelength A' = 
rwindRh, need proper dissipation time > A'/c = r^indRL/c, since the 
current sheets can't expand any faster than the speed of light; alternatively. 
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in a reconnection model, the magnetic field flows into the sheets, to disappear 
in expanding isla nds of hot plasma around 0-lines, with velocity e^Vj^ w ejjc 
( Lvubarskv 12005 ) . with e^j expectecf^ to be on the order of 0.1 to 0.2. In the 
pulsar's center of mass frame, the dissipation time then has the lower limit 

> ^wtnd^L/c; in a reconnection model, Td w r^ind^L/eRC. 




Fig. 16. Magnetic geometry of the inner regions of a striped wind emerging from 
an oblique rotator with a large obliquity i. a) Magnetic structure of the force free 
rotator for i = 60°, from Spit kovsky ( 2006. ). b) One of the two interleaved current 
sheets for the 60 ° rotator, derived from Bogovalov's oblique split monopole model 
jBogovalovl [l999l ) c) The same as b) but for i = 9°, shown for clarity, d) and e) 
Meridional and equatorial cross sections of the striped wind current sheet, for the 
60° rotator, f) A snapshot of a 2D PIC simulation of the equatorial stripes, by 
[Spitkovsk^ (|2002l ) 



An expectation ba sed on kinetic s imulations and experiments on non-relativistic 
reconnection, as in lBirn all (|200lh : relativistic reconnect ion in a pair plasma, the 
case relevant here, has just started to receive attention l|Bessho fc Bhattachariej 
l2005l ). 
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A successful model for the apparent low value of a at the termination shock 
in the Crab Nebula, where the equatorial shock occurs at Rshock ~ 10^i?L, 
requires that the dissipation go to completion in a region where F^md < lO"*'^; 
using the reconnection model reduces this upper limit to F^ind < 10^. In 



MSH15-52, where Rshock - 0.4 - 0.5 pc and Rl = 7825 km, 
a firm upper limit; in a reconnection model, Fwind < 10^'^- 



< 10 



4.7 




J(ion) 0.6 




Fig. 17. Left: Magnetic sandwich geometry of the equatorially concentrated outflow, 
with radial electric current flowing in a extremely weakly magnetized midplane be- 
tween oppositely directed toroidal magnetic field at hig her latitud es. Right: Typical 
magnetic profiles as a function of latitude, from iDel Zanna et all 



It is interesting to note that the full average particle loss rate N± — 
10'*°-5iV4o.5 s-i inferred for the Crab Nebula requires, from energy con- 
servation, that the asymptotic value of the wind 4 velocity is Fwindoa < 
E / N±m±c^ = 10'* '^/7V4o.5 in this system, much less than inferred in ID mod- 
els that neglect th e plasm a required for the radio emission such as that of 
iKennel &: Coronitil ( 19841 ) - the upper limit is achieved if the plasma is cold 
by the time it reaches the termination shock. The acceleration from the in- 
ner wind, launched by the pulsar's magnetic spring to F^ = given in 
pd|) . occurs if the current sheets dissipate and the wind heats, accelerating 
from the resulting interna l pressure gradient. This is a relatively slow process 
(jLvubarskv fc Kirkll2001l ). reaching completion before a fluid element collides 
with the termination shock in the Crab Nebula only if F^ < lO"*'*, N > 10'*° 
s~^. This requires maximal dissipation of the current she ets - in a sheet broad- 
ening model, as originally proposed by (|Coronitilll990h . the sheets must ex- 
pand at a substa ntial fraction of the speed of light in the proper frame of 
the flow. Recently I Arond ( 2007 ) showed that the interaction of the relativistic 
currents in neighboring sheets drives a Weibel-like instability in each sheet, 
with a resulting anomalous resistivity that supports such maximal dissipation 
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in the inner wind of the Crab pulsar, r < 10^ Rr,, a conclusion consistent with 



the maximal dissipation rate model of lKirk fc SkiaeraasenI (|2003l ). However, it 



is safe to say that the transition of the wind from high to low a and from low 
to high four velocity remains not fully understood, and not well constrained 
by direct observation of the winds. 

For a long time, winds have been modeled as having an asymptotic flow 
velocity Fwoo ~ 10^, with various arguments being used, ranging from radia- 
tion modeling of the post-shock flow to the dynamics of high energy particles 
injected by the wind at the shock, with Larmor radii comparable to the ter- 
mination shock's radius, invoked as a dynamical model for the variable wisp 
dynamics near the shock. Such particles, if they exist, have 4 velocities (prob- 
ably gained by acceleration in the wind) much larger than the maximum flow 
4-velocity of the wind inferred from energy conservation. They might occur 
due to acceleration in the current sheet, thus are confined to the equatorial 
sector. 

Some of the energy dissipated in the wind might reappear as an unpulsed 
photon source superposed (because of relativistic beaming) on the pulsar. The 
Crab pulsar exhibits linearly polarized (33%), unpulsed optical emission, with 
intensity ~ 1.2 5% of the main pulse peak intensity with a fixed polarization 
position angle ( Kanbach et aZ.I [ioOSi ) . polarizatio n properties co nsistent with 



emission from the toroidal field in the wind zone ()Barnardlll986f ). If it proves 
possible to interpret such unpulsed flux as radiation from the wind, one might 
obtain significant observational constraints on this difficult problem. 

Wisps as Ion Cyclotron/Magnetosonic Waves 

A non-MHD hypothesis based on kinetic structure in the current sheets' 
plasma (a "beyond MHD" model) does well at reproducing the observed vari- 
ability near the wind termination shcok than MHD schemes. Gallant fc AronsI 



([1994') proposed this current contains a high energy ion beam, accompanied 
by a flux of pairs, with ion energy/particle approaching the total magne- 
tospheric potential energy, and that this ion current carries a large fraction 
of the spindown energy, while still being a minority population by number 
density. 

It is important to note that as far as the dynamics is concerned, high 
energy ions in the equatorial return current, expected to emerge from an 
"acute" pulsar - Z(i7, /x) < 7r/2 - can be replaced by high energy electrons, 
expected to form the return current in an "obtuse" pulsar, Z(17,/x) > 7r/2. 
In both cases, the required acceleration to energy/particle comparable to e'P 
must occur in the win d, perhaps a s particle runaways in the current sheet's 
resistive electric field (lAronsll200i . In the electron case, radiation reaction 



can limit the energy/particle that can be achieved, which makes the acute 
pulsar model slightly preferable, and for the rest of this discussion I confine 
discussion to the ion beam case. 
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Since such ions have Larmor radii comparable to the radius, the compres- 
sions induced in the pair plasma at the ion stream's turning points can appear 
as surface brightenin g spaced with se paration comparable to the spacing of 
the wisps. [Spi tk ovskv fc Arond ( 2004 ) implemented this idea in a time de- 
pendent simulation, showing that in a ID model in a toroidal magnetic field 
winding in one direction in the equatorial flow {i.e., ignoring the reversal of 
the field direction in latitude), the ions' deflection into circular motion in the 
abruptly increased magnetic field at the equatorial shock in the pairs is ion 
cyclotron unstable, with gyrophase bunching forming a compressional limit 
cycle which launches finite amplitude magnetosonic waves in the pairs. This 
launching occurs approximately once per ion Larmor period, which is about 
six months for the parameters determined by comparing the model to the im- 
ages. These waves travel out at speeds « (0.3 — 0.5)c, with the precise velocity 
value depending on the degree of isotropization of the pairs. 

The resulting synthesized surface brightness map looks more than a little 
like the observed waves emitted from the inner X-ray ring i n the Crab Nebula , 
as shown in Figure [HI The observed wave emission period ( Mori et al. 20061 ). 
announced at a conference after the model was developed and published, is in 
good accord with the model's predictions. Application of the model to PSR 
B1509/G320 suggests that "wisp" variability on a time scale of years should 
be found. There is weak, but not ver y convincing, evidence for such variability 
in the partial torus near this pulsar ( DeLaney et al. 20061 ). 



By far the most attractive theoretical feature of the model when it is fit to 
the HST and Chandra movies - fit by eye, there is little use in more elaborate 
fitting procedures, given the ID model's departures from the observationally 
obvious requirement of at least 2D - is the inference that the ion flux required 
to produce the observed surface brightness enhancements is Nion « c#/e, the 
equatorial return current of the force free rotator, in the case of the Crab 
pulsar. Of course, since the magnetic held has largely dissipated in the wind, 
the pairs accompanying the ions must largely neutralize the electric current in 
the ions, but the result is an indication that back at the magnetosphere and 
the stellar surface, some piece of non-force-free electrodynamics does work 
to extract this ion flux so as to maintain the star's charge balance. In turn, 
that suggests i < 90°, although i certainly should be a large fraction of tt/2. 
Another feature of the model is that the fact that the ions in the wind are 
inferred to have Larmor radii comparable to the wind's termination radius 
Rs - this yields ion energy/particle E'ion = e^(mp/me//^e^), and jion ~ 10^'^ 
with the pair multiplicity e valuated in the equator, a value close to the MHD 
wind 4- velocity inferred bv lKennel fc Coronitil (|l984f ). 

The model assumes the underlying acceleration from the neutron star is 
like MHD even in the current sheet, with all the particles - ions and pairs - 
traveling with a single (fluid) 4-velocity until a fluid element encounters the 
shock, even though the flow in question is in the current sheet, where different 
plasma components may have different velocities. Thus Spitkovsky and Arons' 
inference that r^md ~ 10®'^ is based on the assumption that the ions, which 
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carry the electric return current in this model, have the same 4-velocity as 
the underlying and surrounding MHD wind. 




Fig. 18. Ion cyclotron instability model of pulsar wind termination shock variabil- 
ity. Top Left: Chandra snapshot of the X-ray ring and torus in the Crab Nebula. The 
ring is the location of the magnetosonic shock wave in the pairs, whose thickness is 
unresolved. Top Right: Snapshot of a ID hybrid (PIC ions, MHD pairs) simulation of 
the equatorial ion outflow as it encounters the enhanced magnetic field in the pairs. 
The upper half shows the magnetic field with the ion induced compressions. The 
lower half shows the ion phase space {pr /yupatreammionC versus r/Rshock, Pr = ra- 
dial momentum of an ion), with the gyrophase bunch formed as the ions encounter 
enhanced B and acquire a reduced Larmor radius. The rotating, reforming bunch 
(a limit cycle) emits a train of compressional magnetosonic waves that propagate 
to larger radius. Bottom: Snapshot of the synthetic surface brightness map, assum- 
ing the pairs emit synchrotron radiatio n in the compressed magne tic field with no 
nonthermal acceleration included. From ISpitkovskv fc Arona (|2004l l. 
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The upper limit F^md < lO'*"'^ clearly is inconsistent with a single fluid 
velocity for ions and pairs. The ID model averages over the whole equatorial 
sector that feels the equatorial belt shock, shown in Figure[T31 thus mixes the 
ion flux with the whole flux of pairs feeding the torus, ~ iQ'^* ''^ pairs/sec for the 
optical, X- and 7-ray emission from the Crab. It yields an average value of a 
within a factor of 3 of the average value of a inferred from the MHD nebular 
models. However, the assumption Fwind = lion clearly violates the upper 
limit on Fwind, a limit which comes from assuming the whole flux of pairs 
feeds the equatorial toruJ^. an ass umption supported by the approximate 
correspondence of the radio wisps (iBietenholz et al. with the optical 

features. 

The sandwich magnetic field m odel, a fundamentally 2D c onstruct, sug- 
gests substantial alterations of the Spitkovskv fc Arong ( 2004 ) scheme. The 
magnetic field inside the current sheet is weak compared to that in the MHD 
flow outside the sheet. That weakened B alters the character of the ion orbits 
from simple magnetic reflections to partial deflections from radial flow, thus 
altering the momentum transfer to the pairs and therefore the compressions. 
If most of the ions flow in an essentially unmagnetized region, 7io„ is no longer 
coupled by the magnetic fi eld to the flow 4- velocity of the pairs, thus allowing 

lion -^uimd 

(lArondl2007t ). 

Taking such improvements of the model into account is needed before one 
can realistically assess the model's consequences for observations, including 
possible hadronic TeV gamma ray and neutrino emission from the nebulae, 
a definite prediction of the model at some flux level. At this writing, mod- 
els of the TeV SEP based on the ID dynamical model ( Amato et a\ 2003t 
Bednarek fc Bartosik |2003^ show that the hadronic gamma rays at e < 20 



TeV are masked by inverse Com pton emission; recent HESS observations of 
the Crab ( Aharonian et g/.ll2006[ ) suggest t hat an interes t ing co nstraint on ions 
in the outflow may be attainable. Recentlv lHorns et al. ( 20061 ) suggested that 
the TeV emission observed from the Vela-X PWN might be due to hadronic 
emission from the p-p interaction. Simple evolutionary models (Bucciantini 
and Arons, in preparation) suggest that such an interpretation is supported 
by the ion current sandwich model, when proper account is taken of the ex- 
pansion and compression history of this older (age ^ lO'*'^ years) PWN. 

Whatever the ultimate fate of models of this sort, they suggest the use- 
fulness of stepping beyond MHD (which does not by itself constrain the wind 
velocity or density) in modeling the observations with the goal of extracting 
the plasma content and constraining just what does come out of pulsars, and 
other compact objects. 



lArond (|l998al ) suggested the large flux of radio emitting particles might be in 
higher latitude flow, a possibility I now disfavor, both because continuity of the 
Crab Nebula's spectral energy distribution militates against the accident that 
would be required, if the particle spectra were formed from such diflferent flow 
components, and because the modern MHD models have no such distinction 
between low and high latitude particle fluxes. 
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See lKirk. Lvubarskv fc Petril (|2007f) for a parallel review of this subject. 



Shock Acceleration 



Starting with the seminal work of iRees fc Gund (|l974l ). the conversion of 



flow energy to the nonthermal particles emitting nebular synchrotron and 
inverse Compton radiation has been attributed to some form of "shock ac- 
celeration". Most commonly, Diffusive Fermi Acceleration (DFA) has been 
invoked, even though in relativistic shocks it faces a number of substantial 
difficulties, especially when the magnetic field is transverse to the flow. Mech- 
anisms h ave not been apparent to su pply the very large amplitude turbulence 
required ( Niemic fc Ostrowski 20061 ). which must extend to large distances 



('^ Rshock) both up and downstream of the shock so as to have a large "opti- 
cal depth" for scattering of the largest Larmor radius particles both up- and 
down-stream. The spectrum of test particles accelerated depends sensitively 
upon the isotropy (or lack thereof) of the scattering process. However, in the 
case of isotropic scatterin g in the fluid fram e, the test par ticle spectra are en- 



case ot isotropic scatterin g m tne iluiQ tram e, tne test par ticle spectra are en- 
couraging - Monte Carlo (jKirk et a/.ll200o[ ) and analytic ( Keshet fc WaxmanI 



|2I3G5) calculations with assumed scattering rules and infinite optical depth 
for particle scattering up and downstream yield an accelerated particle spec- 
trum in the downstream medium N[E) oc i?^^"/^, which is almost exactly 
that inferred by modeling the synchrotron emission in a ID post-shock flow 
in the Crab Nebula ([Kennel fc Coronitil Il984h - this simplified flow model 
should be a not unreasonable approximation to flow right outside the equato- 
rial belt shock shown in Figure 1141 The efficiency depends entirely on what is 
assumed for the particle injection rate into the process, and the acceleration 
rate depends entirely on the assumed turbulence amplitude that goes into the 
sca ttering law adopted. 

Hoshino et "oil (|l992h suggested an alternate process, especially well tuned 
to the mixture of heavy ions and pairs injected in the equator with the mag- 
netic field tra nsverse to the flow^i. Using ID PIC simulations, they, and, 



more recently, Amato fc Arond (j2006l ) with higher mass ratio and resolution. 



showed that high harmonics of the ion cyclotron waves generated by the ions 
as they pass through the shock in the pairs can be resonantly absorbed by the 
positrons and electrons, which are non-thermally heated, yielding power law 
downstream distributions with a spectral slope that depends on the ratio of 
the upstream ion energy density to that of the pairs. 

The nonthermal part of the spectra shown in Figure [19] extends from 
the pairs' flow energy /particle all the way up to the ions' energy/particle. If 



As in the high energy beam, cyclotron instability interpretation of the wisps, an 
electron current accelerated in the equatorial current sheet of an "obtuse" pulsar 
can replace the ion beam accelerated in the current sheet of an "acute" pulsar 
without altering the conclusions - at equal relativistic energy/particle, the only 
difference is the sense of gyration with respect to the unknown vector direction 
of the magnetic field. 
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all the species have the same upstream flow velocity, the resulting spectra 
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nicely span the range requ ired for optical, X-ray and 7-ray emission from the 
Crab (|Hoshino et aLlll992l ). and for X-ray and 7-ray emission from G320 
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Fig. 19. Downstream particle spectra of a PIC simulation of a shock with upstream 
magnetic field transverse to the flow in an electron-positron-proton plasma with 
mass ratio mp/m± = 100 and upstream flow Lorentz factor of all species 71 — 40. 
The upstream Lorentz factor sets the scale; the results are otherwise independent of 
the specific value of 71, so long as it is larger than 2 or 3. /{y) are the distribution 
functions, with J /(7)7rf7 = species density. The ratio of the proton number density 
to the number density of electrons plus positrons was 1/19, therefore the energy 
density in upstream protons was 2.1 times the energy density in pairs. The dashed 
curves in the panels showing the electron and positron spectra are the best fit 
Maxwellians, while the temperatures stated in the figure (which are in units of m±c?) 
rows are inferred from fitting a Maxwellian with a power law tail to the numerical 
data. The slope of the power suprathermal component is a = 3.2, corresponding to 
an energy space spectrum dN/dj oc i"'""^^ oc 7~^'^, quite close to what is required 
in modeling the optical and X-ray synchrotron emission from the Crab Nebula and 
other young pulsar wind nebulae. The power law spectral index is a strong function 
of the upstream energy density ratio, thus making the downstream nonthermality of 
the pairs a strong function of the upstream composition. For the details, including 
the p ower law spectra index as a function of density ratio, see lAmato fc Arona 
ll2006t) . 
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( Gaensler et al. 2002 ). As demon strated by Amatq fc Arond ( 2006^ . however, 
when apphed to the ID model of lGallant fc Arond (|l994r ) of the Crab Pulsar 
wind's termination shock, this mechanism has trouble providing an accelerated 
spectrum of pairs in accord with the observation. 

The ion flux by number is fixed at Goldreich- Julian value (it can hardly be 
anything else). If 7ion = ^mmd, the upstream energy density ratio {U±/Uion)i = 

) ~ 10^ leads to the pairs' particle energy distribu- 
tion hardly differing from a relativistic Maxwellian, the do wnstream distribu- 

tion for a relativist ic transverse shock in a pure pair p lasma (jLangdon. Arons fc Max 
Ii988; G allant et al . 1992; Spit kovskv fc Arondl2007n . not at all in accord with 
the observations. 

If acceleration occurs near the sandwich midplane, the obstacles to accel- 
eration by either mechanism may be reduced. The cyclotron mechanism bene- 
fitsfrom 7io„ > r^atnd- With pair multiplicity k± = 2n±/nGj > lO*" (required 
to supply the radio emission of the Crab Nebula) and therefore < 
(since most of the energy flux is carried by the pairs) while 7ion ~ 10^'^, now 
{U±/Uion)i ^ 3, which leads to a downstream particle spectrum possibly as 
flat as the E~^^^ radio emitting particles. DFA might benefit from a weaker 
transverse magnetic field - for criocai < 10^"^ within the current sheet, mag- 
netized shocks in pair plasmas become indistinguishable from sho cks formed 
in a flow with no magnetic field at all ( Spitkovskv fc Arond 2007 ) . Figure [20l 
shows a snapshot of a 2D PIC simulation of a shock in a B = 0, plasma, 
exhibiting a test particle gaining energy as it scatters in the magnetic turbu- 
lence in and behind the shock front, which forms due to the Weibel instability 
driven by the upstream flow penetrating into the heated downstream medium. 

The downstream pair spectra found in this simulation are shown in Figure 
[2T1 which exhibit the formation of a suprathermal particle spectrum through 
scattering in the turbulence in and near the shock front. To go to energies 
much higher requires turbulence that persists to much greater depths in the 
downstream than are studied in the simulations, with amplitudes that do not 
decay. Phase mixing betwee n particles and fields may cause th e decay of the 
downstream magnetic fields ( Chang. Spitkovskv fc Aroii3l2007 ). in completely 
unmagnetized shocks. Nevertheless, these results, which show how nonthermal 
particles can be created out of the thermal pool in a very weakly magnetized 
shock (upstream a < 10"'^, as might be characteristic of the central regions 
of the equatorial current sheet), there to act as seeds for DFA. Identifying the 
necessary scattering turbulence remains a challenge. 

These speculative possibilities suggest a solution to the vexing question 
of how the peculiar spectrum of particles injected into the Crab Nebula and 
other PWNe might be formed. These systems all have very flat power law 
distributions of particles N{E) (x E~p, 1 < p < 1.5 giving rise to their radio 
synchrotron emission, while the inferred injection spectra of particles required 
for the harder photon emissions (optical when seen. X-ray) have much steeper 
spectra, p — 2.2 — 2.4. The continuity observed between the radio and the 
harder photon spectra suggests the shock injects a broken power law spec- 




Fig. 20. Shock structure and test particle gaining energy at the sho ck front, from 
a PIC simulation of an unmagnetized shock by ISptikovskvl (|2007h . Lengths are 
measured in units of the upstream skin depth. Bottom panel: Density as a function of 
position, exhibiting the factor of 3 jump (properties of this shock in a 2D plasma are 
measured in the downstream frame). Middle panel: {B^y^^, showing the upstream 
magnetic filaments characteristic of the Weibel instability in the linear regime, the 
scrambled magnetic structures formed when the currents reach the Alfven critical 
curre nt, magnetic trapping disrupts the flow and the shock transition forms (|Katcil 
I2OO5I ). Note the transition of the magnetic structure to spatially intermittent (widely 
separated) islands in the downstream. The orbit of a test particle is superimposed. 
Because of the finite thickness of the strongly turbulent scattering layer, particles 
escape downstream after a finite number of scatterings. The sideways dimension is 
periodic, thus a particle leaving the box at the top reappears at the lower bottom. 
Top panel: Energy of the test particle, which started with the upstream flow energy 
(7 — 15) and increased its energy by a factor of 10 before escaping. 



trum. Cyclotron resonant acceleration, driven by an ion flow in the current 
sheet with jion ^ , might be responsible for the very flat spectrum of radio 
and infrared emitting electrons, whose upper cutoff is determined by energy 
conservation to be not large compared to F^md- DFA might be responsible 
for continued acceleration to optical. X-ray and gamma ray synchrotron emit - 
ting energies. This combination was first suggested bv lGallant et MI (l2002h . 



In their formulation, in which 7io„ = F^ind, the ions had negligible energy 
density compared to the much denser pairs and therefore could not act as the 
desired accelerator, while in the version suggested here, the ions having a 4 
velocity much larger than that of the pairs allows them to act as an effec- 
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Fig. 21. Downstream particle distribution function in the large PIC simulation of 
a pair shock with no systematic magnetic field in the upstream medium used to 
create Figure [20l The lowest dashed is a Maxwellian distribution fit. The solid curve 
is the actual particle distribution. Clearly, scattering in and around the shock front 
produces suprathermal particles (in this example, these are 5% by number and 20% 
by energy). 

tive agent in creating the nonthermal radio emitting particles. A quantitative 
expression of these thoughts is under investigation. 

The scheme outlined above assumes the stripes in the wind do in fact 
disappear well upstream of the shock, as is suggested by the rapid dissipation 
results of|Arons (2007). Lyubarskv (2003, 2005), starting from the slow stripe 
dissipation model of iLvubarskv fc Kirki (j200l[ ). showed, using ID and 2.5D 
kinetic simulations in plasma, that the striped field can annihilate in a 
broad "shock" region with strong plasma heating - thus, the effect is as if a 
was low in the upstream medium, even when it was not in fact small. He also 
presented evidence that reconnection at individual sheets might create a flat 
particle distribution N cx E~^'^ which might be of use in understanding the 
radio emission from PWNe. However, a full 3D treatment is needed in order to 
properly evaluate the nonthermal particle acceleration, since the periodicity 
of the simulation in the electric field direction can lead to artificial results for 
particle acceleration. 

A clear theoretical picture has yet to emerge, never mind models with 
readily checkable observational predictions. And all of these models and the 
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observations underlying them tell us that a fully quantitative account of the 
pair plasma supply is still lacking. There will be progress on the theoretical 
front in the next few years. The much anticipated launch of GLAST may add a 
new observational handle on these problems, since the highest energy radiating 
particles have rapid radiation losses, therefore leading to interesting time series 
in the gamma ray emission as the particles simultaneously accelerate and 
radiate in the time variable termination shock region. 

4 Conclusion: Pulsar Problems and Prospects 

Pulsar physics has made substantial progress in the last decade. The emer- 
gence of the MHD model of rclativistic plasma flow in Pulsar Wind Nebulae 
has given a plausible account of the plume-torus (a.k.a. jet-torus) structures 
observed in these systems, thus reconciling the external response of the sur- 
rounding world to the well known difhculties in forming a jet in the rclativistic 
winds themselves. These models also explain the lack of shock excitation of 
the nebulae at high latitude as a result of the shock's lack of spherical sym- 
metry, itself a consequence of the anisotropy of the energy outflow found in 
the energy flux emerging from the magnetosphere. 

Application of force-free electrodynamics and rclativistic MHD to the 
winds emerging from magnetospheres with dipole magnetic fields has led to 
the first theory of the oblique rotator's energy loss that incorporates qualita- 
tive changes from the vacuum theory imposed by electric current flow, leading 
to the delightfully simple formula for the spindown energy loss given in ex- 
pression ([6]). For the model of the magnetopshere, the most important result 
is that the polar electric current distribution is close to that of the monopo- 
lar magnetopshere, reflecting the asymptotically monopolar poloidal magnetic 
field beyond the light cylinder. 

The role of reconnection in the transfer of open to closed magnetic flux 
(and back again, since this is an unsteady process) has begun to be assessed, 
and is full of promise as a path to a physical theory of the boundary layer 
between the closed and open magnetosphere, where existing gamma ray obser- 
vations and gap models suggest the most prominent photon emissions from 
pulsars occur. Prominent issues waiting assessment include reconciling the 
creation of pairs with the monopolar current distribution in expression ([9|) 
and Figure [HI with various solutions being on the table, awaiting surgery - 
these range from rapid local current fluctuations (averaging to the force-free 
current) to manipulation of the polar cap electrostatics by return currents 
on the open flux tube boundaries. These models all have consequences for 
long standing issues such as the origin of torque fluctuations, radio subpulse 
phase randomness and drifting, possibly for the origin of the departures of 
the braking index from its canonical value of 3, and for the origin of the large 
particle fluxes inferred from nebular radio emission. The modeling will be ob- 
servationally illuminated by the results of the upcoming GLAST gamma ray 
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mission. If high sensitivity X-ray astronomy has a future, observations with 
the abiUty to inspect variabihty in the X-rays, both nonthermal and thermal 
from polar caps heated by magnetospheric currents would be invaluable. 

The longstanding problem of the origin of the weak magnetic fields inferred 
downstream of pulsar winds' termination shocks in the young nebulae is still 
an outstanding question, with dissipation of the magnetic stripes in the wind 
being the prime suspect. Whether this occurs in the wind far upstream form 
the shock, or in the shock itself, is an open question under active investigation. 

Finally, the basic physics of relativistic shock waves is receiving signifi- 
cant attention, which opens the prospect of having a testable theory of the 
conversion of flow energy to non-thermal particle spectra in these relativistic 
systems within the next few years. Significant issues that will receive atten- 
tion include the mixture (if any) of Diffusive Fermi Acceleration and other 
shock related processes, the possible role of protons and other heavy ions as 
well as pairs in the acceleration physics (and testing for these ions' presence 
through VHE gamma ray and neutrino observations), and an understanding 
of how the nebular radio emitting electrons can be accelerated by the curved 
termination shock. 

Exciting times lie ahead! 
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